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FOREWORD 

This report presents the results of work performed by 
Lockheed Missiles & Space Company, Inc., Huntsville Research 
& Engineering Center, under Contract NAS8-28899 for the NASA- 
George C. Marshall Space Flight Center, Marshall Space Flight 
Center, Alabama. The performance period covered by this final 
report is from 30 June through 31 December 1972. Mr. Charles 
S. Cornelius of the Control Mechanisms Branch of the Astrionics 
Laboratory's Guidance and Controls Division is the Technical 
Monitor for this contract. 
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Section 1 

INTRODUCTION AND SUMMARY 

Space vehicles and missiles which have active guidance systems require 
a means of implementing control of the vehicle. This control function may be 
required over large ranges of velocity, acceleration, altitude and vehicle 
orientation (angles of attack, yaw and roll). There are four basic means for 
providing this control: (1) inertial systems such as those used on satellites 

in orbit; (2) reaction control systems (RCS) which utilize small thrusters; 

(3) aerodynamic control surfaces used in the sensible atmosphere; and (4) 
thrust vector control systems (TVC) which change the direction of the engine 
thrust vector. The means of control used on a vehicle is dictated by the flight 
regions in which it operates. 

Future space vehicles will fly ascent and reentry trajectories which 
pass from the sensible atmosphere to near vacuum environment and 
return. These vehicles will therefore use a combination of aerodynamic sur- 
face controls, TVC or RCS. It is anticipated that the TVC system will be 
relied upon to provide launch vehicle stability and control over much of the 
flight trajectory. One type of TVC is secondary injection TVC. 

The characteristics and performance of secondary injection thrust vector 
control (SITVC) systems are highly dependent upon the hardware character- 
istics used to implement the system, the injector geometry, nozzle injection 
locations, injection attitude, main nozzle characteristics, main propellants, 
and secondary injection fluid. Scale model tests of secondary injection 
systems can provide useful data for the evaluation of SITVC. However, 
depending upon the systems being modeled, significant variations in the degree 
of simulation can be encountered unless attention is given to the gasdynamic 
and kinetic scaling of the interacting flows. 
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The objective of this effort is to evaluate a scale model of a low-pressure 
rocket engine which is to be used for secondary injection studies. The following 
specific tasks were pursued under this contract: 

1. Assess the test conditions which are required to achieve full-scale 
simulations. For scale model tests where full-scale simulation is not achieved, 
perform calculations to obtain estimates of the effects of non- simulation on 
SITVC performance. 

2. Recommend fluids to be used for both primary and secondary flows, 
and suggest test operating conditions such as chamber pressures, injection 
pressure, injection velocity and other relevant parameters. 

3. Recommend possible modifications to be made to the scale model 
and its test facility to achieve the highest possible degree of simulation. 

This report presents a discussion of the theoretical and empirical 
scaling laws which must be observed to apply scale model test data to full- 
scale systems and describes a technique by which the side forces due to 
secondary injection can be analytically estimated. 

Based on method-of-characteristics calculation for a typical full-scale 
nozzle and its scale model, recommendations are given with regard to scale 
model nozzle geometric and test fluid parameters which must be observed in 
order to suitably simulate the SITVC system effects. Also included is a 
critical evaluation of the scale model test facility being considered by MSFC. 


2 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D306359 


Section 2 

TECHNICAL DISCUSSION 

2.1 CONCEPTUAL BACKGROUND OF THRUST VECTOR CONTROL 

The concept of thrust vector control by secondary injection was con- 
ceived in 1949 by A. E. Wetherbee, Jr. (U. S. Patent 2,943,821), Ref. 1. Much 
research and development effort has been conducted in this area since 1952, 
and the concepts have been extended to cover supersonic mixing and combus- 
tion. Most of the investigations have been conducted for the case of a two- 
dimensional slot or round jet on the surface of a two-dimensional flat plate 
expanding into a uniform approach flow (Refs. 2 through 6). Ambient temper- 
ature air or inert gases were used in most of these studies. Although the 
results of these studies have been extremely important in the development 
of the TVC concept, they are limited to basic two-dimensional flow applica- 
tions. Future applications necessitate that the TVC state of the art be 
extended to include an understanding of the jet interacting with complex, 
three-dimensional approach flow. 

2.2 SCALING PARAMETERS FOR TVC STUDIES 

The task of scaling the jet interaction phenomena has been practiced by 
many investigators as evidenced by the large number of reports on experi- 
mental studies. Although much knowledge has been gained by these works and 
the techniques for scaling the jet interaction phenomena are well established, 
it is often difficult to duplicate the scaling parameters from a practical aspect. 
Therefore it is of the utmost importance to examine the scaling techniques that 
have been used in the past and evaluate them in terms of future applications. 

As stated previously, the jet interaction phenomena as defined for this 
program are composed of the complex interactions of a sonic jet expanding 
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from the nozzle wall with the supersonic main flow surrounding the secondary- 
jet. Scaling these phenomena for a prototype engine involves scaling the engine 
nozzle and secondary jet nozzle geometry, considering the properties of jet 
working fluid, and gas dynamic characteristics of jet and main nozzle flows. 

Geometric scaling of the prototype engine is an accepted practice for 
model tests and is of particular importance in terms of the jet interaction 
effects. The relative size, location and orientation of the jet port on the proto- 
type must be maintained on the model. This is to ensure that scaled jet plume 
interaction with adjacent surfaces will be possible, assuming that the correct 
gas dynamic simulation of the jet has been achieved. Therefore, in a para- 
metric study of the jet interaction effects, the sensitivity to jet size, location 
and orientation should be investigated. 

In gasdynamic scaling of rocket engine flows, matching of the Reynolds 
number is an accepted criterion. Recalling the definition of Reynolds number 
as puL/p, it can be shown that for a geometrically scaled model, using the 
full-scale gas, the Reynolds number of the flow in a nozzle can be matched 
only if the density, p, is increased via a chamber pressure (P^) increase as 
the model scale is decreased. (For a one-tenth scale model, a chamber pres- 
sure equal to ten times full scale is required to match full-scale Reynolds 
number.) If the required Reynolds matching can be achieved the resulting test 
data should give full-scale results (presuming full-scale operating fluids are 
used). 

Relating the above information to SITVC studies, it becomes immediately 
apparent that the Reynolds numbers of both the injector flow and of the flow in 
the primary nozzle should match full-scale values. It is determined rather 
quickly, however, that the above conditions are not practiced for most test pro- 
grams. The magnitude of the pressure levels required to obtain full-scale 
Reynolds number in scale model TVC programs are, in general, beyond the 
practicalities of existing test facilities. It is also quickly determined that the 
use of the full-scale gases may be impossible or impractical due to handling 
problems or the difficulty of duplicating the combustion products at a reason- 
able cost. 

4 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D3 06359 


Accepting the fact that full-scale Reynolds numbers of the primary and 
injector nozzle flows cannot be matched, the next most desirable condition is 
that which would minimize the effects of Reynolds number. This condition is 
generally accepted to exist when the local Reynolds number is of sufficient 
magnitude such that the flow adjacent to the body is turbulent; i.e., the local 
Reynolds number is greater than or equal to 10^. 

Methods for determining boundary layer transition is a separate study 
and will not be discussed here; however, previous analysis of large rocket 
engines has confirmed that the flow in the region of the TVC injection ports 
is fully turbulent. The following discussions of the problem of simulating the 
full-scale system will be predicated on the assumption of the model operating 
in the turbulent flow regime. 

From the preceding discussion it becomes immediately apparent that 
"direct" scaling cannot be achieved in most TVC test programs. A "simula- 
tion" technique must therefore be employed where the effects of Reynolds 
number are considered to be second order and negligible. This results in the 
dilemma faced in most TVC model testing — what parameters must be con- 
sidered to "simulate" the full-scale conditions and how does one accomplish 
the simulation ? 

In general, it is extremely difficult to simultaneously simulate (theoret- 
ically or experimentally) both the inviscid and viscous aspects of the gas dy- 
namics. In some instances simulation of only the inviscid characteristics has 
been all that is required to obtain meaningful results. It is anticipated that 
meaningful simulation of the flow phenomena to obtain TVC pressure informa- 
tion will require consideration of both inviscid and viscous effects. The 
simulation of these effects is discussed in the paragraphs that follow. 

The complexity of the TVC flow problem has obstructed the development 
of a general analytical solution and increased the need for developing experi- 
mental techniques. To obtain meaningful experimental results, the full-scale 
flow field must be simulated. The parameters required for similitude must, 
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therefore, be identified and substantiated through experiment. Two phe- 
nomena associated with TVC performance which, should be simulated are 
the penetration depth and the plume shape resulting from secondary injection 
into the primary flow. 


The penetration of the injector flow into the primary stream is essen- 
tially one case of the expansion of a sonic or supersonic jet into a nonquiescent 
stream. It would therefore seem logical that simulation parameters evolved 
for similar problems might also be applicable to the TVC problem. Goethert 
and Barnes (Ref. 7) addressed the problem of simulation of base flow phe- 
nomena with a cold gas. They showed that for simulation of a supersonic jet 
expanding into a parallel supersonic flow, the following inviscid parameters 
should be satisfied: 


Geometric scaling of the nozzle 
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Satisfactorily matching these parameters results in the matching of the 
shape of the full-scale inviscid plume boundary. In order to match not only 
the inviscid plume boundary, but also the detailed inviscid flow structure 
(which is necessary for achieving correct inviscid force ratios) then it is also 
mandatory to have the correct values of y and y . For the case of the TVC 
studies, matching these parameters should result in the correct inviscid re- 
lationship between injector and primary flows. 


Since the basic nature of the interaction of a jet with an external flow in- 
volves viscous effects and the mixing process, the relations which describe these 
phenomena were reviewed to establish conditions for viscous similitude between 
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the full-scale and model flows. To determine a "viscous mixing" similarity 
parameter, empirical eddy viscosity models of several types were applied to 
the momentum equation for steady, compressible flow (Ref. 8). The models 
used and the resultant similarity parameters are summarized below: 
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Based on the results of Ref. 8 it was concluded that matching of the ratio of 
the "pu" product for injectant and primary flows between prototype and model 
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is an apparent similarity parameter for -TVC flow simulation. A systematic 
parametric test program will be required to establish the pu ratio as a 
similarity parameter. 


The "pu" ratio parameter can be rewritten in the following form using 
the equation of state and the definition of sonic velocity: 
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An examination of this expression reveals that the "pu" ratio parameter im- 
poses an additional requirement on the ratio of the sonic velocity of the in- 
jectant and primary gases. Rearranging the above expression and assuming 
that the "yPM" ratio for prototype and model are equal (i.e., 
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Therefore, in summary, an examination of the requirements for gas- 
dynamic similitude of full-scale and model flows has resulted in the following 
similarity conditions: 


• Geometry — full-scale primary and injector nozzle contour 
geometry properly scaled. 

• Gasdynamic scaling 
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These parameters can be readily satisfied if gamma (y .) of the model jet is 

J 

equal to y. of the prototype. If y s are not equal, then adjustments must be 
J J 

made in nozzle geometry to change the model area ratio and ultimately the 
jet chamber pressure. The parameters P^/p^ and y PM^)^/y PM have 
been used in Refs. 10 and 11, respectively, to correlate analytical and experi- 
mental jet interaction data. Inherent in the above parameters is the require- 
ment that the primary nozzle flow conditions be the same for model and 
prototype. 


By appropriately combining the above equations Zukoski and Spaid 
(Ref. 12) have determined a scaling parameter convenient for TVC applications. 
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This scaling parameter can be stated as 



prototype model 

The appearance of total temperature, T , and molecular weight, W, intro- 
duces additional complications. In addition to satisfying the above relation- 
ship, the jet nozzle must be geometrically scaled, and the ratio of jet to free 
stream total pressure (P q /p q ) must be maintained to ensure no scaling 
errors. J 00 

Correlation of interference pressure data for different pressure ratios 
and nozzle diameters for an underexpanded jet in a supersonic flow was 
obtained by using a length parameter, h, in Ref. 13. This parameter denotes 
the distance from the nozzle exit to the Mach disk or to the first intersection 
of the "diamond" shock pattern when this configuration exists. The parameter 
"h" is defined by the empirical equation: 



For correlation purposes d is replaced by d t> jet nozzle throat diameter. 

For conditions where the exit Mach number of the model and prototype jet 
nozzle are the same, the above equation can be applied to scale the jet inter- 
action phenomena. Examination and manipulation of these correlation/scaling 
parameters can show that they are inherently included in the gasdynamic 
scaling parameters summarized on page 9. 
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Therefore, the requirements which must be met in order to accurately 
simulate the secondary injection phenomena for a non-reacting system can be 
summarized as follows: 
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The sensitivity of the scaling of the jet interaction phenomena to varia- 
tions in the basic gasdynamic parameters is not known. The applicability of 
scaling parameters evolved from simple two-dimensional models is also in 
question. Therefore, without any obvious optimum choice, scaling of jet inter- 
action phenomena for this study will be based on the similarity parameters 
which have been used to simulate the plume of an underexpanded jet exhausting 
to a supersonic flow. In addition, investigators have reported dependency of 
jet interference on the injectant molecular weight and total temperature (Refs. 
14 and 15). Therefore, the sensitivity of the jet interaction phenomena to vari- 
ations in the following gasdynamic parameters should be investigated: 


• Jet exit to engine flow static pressure ratio, P./p 
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• Jet exit gamma, y. 

• Jet exit molecular weight, W. 

• Jet total temperature, T 

o . 

J 

• Engine flow approach conditions, M , P 

oo oo 

In the process of conducting the tests in which the above gasdynamic 
parameters are used as independent variables, the effects of the scaling 
relations will also be obtained. By selecting the proper test fluids, each of 
the above gasdynamic parameters can be varied independently during the 
course of a parametric study. 

2.3 SIDE-FORCE PREDICTION TECHNIQUES 

The phenomenon of a secondary jet expanding perpendicularly into a 
supersonic mainstream has been described and documented by various investi- 
gators, for example Refs. 12, 13, 16, 17 and 18. The generally accepted major 
features of this phenomenon are schematically defined in Fig. 1 of the Appendix. 
The turning of the nozzle flow due to interaction with the secondary jet sub- 
stantially raises the pressure (P 2 ) on the surface of the nozzle wall upstream 
and immediately around the jet port(s). The integrated effect of these pressure 
increases over the separated flow region augments the pure reaction force of 
the jet. The ratio of the sum of this integrated pressure force (F^) and the 
TVC motor thrust normal to the engine surface at the point of injection to the 
jet motor vacuum thrust (F^) defines the TVC force amplification factor, AF, 
as follows: 


AF 


F /F 
y' x 

m /m 
s P 
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At realistic conditions with the main engine firing, the amplification factor 
can increase up to approximately four, depending upon the approach flow prop- 
erties, secondary jet gas properties and nozzle geometry. 

The jet interaction phenomenon has been mathematically modeled by 
many investigators. A summary of some of the more pertinent models is 
shown in Table 1. These models are, in general, limited to two-dimensional 
or modified two-dimensional flow situations and are not applicable to problems 
associated with the flow in contoured rocket nozzles which is three-dimensional. 

The development of a three-dimensional slot jet mathematical model 
which addresses some of the problems associated with TVC application was 
accomplished in a Lockheed Independent Development Study during 1969 and 
1970. This model, which is capable of estimating supersonic boundary layer 
jet interaction control force characteristics, has been coded for application 
to digital computers and is available for engineering studies. A brief discus- 
sion of this model is presented in the following paragraphs. A more detailed 
description of the development of this model is included in the Appendix of this 
document. 

The three-dimensional gas interaction mathematical model, Fig. 4 of 
the Appendix, is based upon a momentum balance. It consists of equating the 
drag on the slot plume (equivalent) to the change in axial momentum of the 
transverse jet. The basic assumptions used in formulating this model are 
summarized in Table 2. The mean reattachment pressure on the upstream 
face of the injected plume and the base pressure on the downstream side, in 
the present computer program, are computed using Spaid and Zukoski' s 
empirical constants (Ref. 4). However, if desired, the downstream empirical 
constant can be computed by using the method proposed by Maurer (Ref. 5) to 
handle the effects of slot inclination angle while the upstream empirical con- 
stant may be related to the mixing theory proposed by Korst (Ref. 18). Flow 
which turns around the end of the slot and proceeds downstream is assumed 
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to have a boundary layer separation pressure value over a quarter circle 
❖ 

area with a radius equal to the computed separation distance at the end of 
the slot span. The end effect significantly increases the boundary layer 
separation force produced by the jet interaction phenomena. The magnitude 
of the ratio, boundary layer separation force to jet reaction force, is a func- 
tion of: local freestream Mach number, slot aspect ratio, injected gas total 
to freestream static pressure ratio, nozzle geometry, the direction of nozzle 
inclination, and the gas properties associated with the two streams. The 
Lockheed program will compute the penetration height of the injected gas and 
the point of freestream flow separation in front of the slot as a result of three- 
dimensional flow over and around a transverse jet injected into a freestream 
flow field. Although the above is applicable for a flat plate, the incremental 
force terms can be warped to account for body curvature within the existing 
computer program. 

The INPUT and OUTPUT parameters available in the Lockheed program 
are presented in Table 3. The sensitivity of the amplification factors to pro- 
gram input parameters is discussed in Table 4. 

2.4 GASEOUS SECONDARY INJECTION SIDE FORCE CALCULATIONS 

The procedure presented in this section for calculating thrust vector 
control performance was developed by personnel of the Aerotherm Corporation 
(Ref. 19) as a part of Air Force Programs AF04(6l l)-9075 and AF04(6l 1 )-9960. 
The technique appears to be applicable to predict both liquid and solid motor 
TVC system performance as shown by the data comparisons in Dahm' s report 
(Ref. 19) 


A blast wave solution has also been programmed and checked out on the 
computer; however, further comparisons with test data are needed to verify 
its accuracy. For this reason and since the above theory applies to slot 
injection the approach of Section 2.4 will be used. 
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The basic approach for developing a suitable calculational technique is 
first to determine the nature of the components which make up the TVC side 
force, then to arrive at equations which describe these forces, and finally to 
ascertain how these forces make theirselves felt in the form of vehicle turning 
moment. The subsequent paragraphs describe each of these items in detail. 

The injection of a non- reacting gas into the supersonic exhaust flow of 
a rocket nozzle is a means to generate lateral force on the nozzle with a re- 
sulting turning moment to provide pitch control for the vehicle. The lateral 
force is a result of the contributions of the momentum thrust of the injection 
flow and the interaction forces due to the separation of the primary nozzle 
flow by the secondary flow. Figure 1 presents a component diagram of these 
TVC forces (single port injection). 

The equations to calculate the magnitude of the TVC forces are in part 
based on the following assumptions: 


a. The state of the primary gas may be calculated using real or 
ideal gas thermochemistry. These calculations will be made 
using a recognized technique e.g., the method of characteristics. 

b. Flow within the secondary injection port is assumed to be 
described adequately by a one-dimensional analysis. The 
secondary nozzle is also assumed to be circular at its exit 
with the flow considered to be adiabatic. 


The forces due to the momentum of the secondary jet are composed of 

forces acting parallel to the secondary nozzle centerline (F ) and normal to 

the secondary nozzle centerline (F ). 

1 mn 

The basic one -dimensional momentum equation describing the thrust 
generated parallel to the secondary nozzle axis is: 


F 

m 



(p' A' u' 2 + p A ) 

^s s s *s s ideal 



A 

s 


( 1 ) 


15 

LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D3063 59 


where the primed parameters refer to the one -dimensional value of these 

parameters at the secondary nozzle exit and p is the secondary nozzle 

s 

thrust efficiency. The value of r) accounts for the thrust losses due to 

n 

s 

surface shear forces, non one -dimensionality and real gas effects. Dividing 

$ 

all terms in Eq. (1) by the actual "corrected" mass flow rate, m C one 

s s 

arrives at a form of the equation which is more convenient for TYC calculations 


m 
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where C^. comes from the following equation 
v 
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and 


and 


r = 




y + 1 
2(7- 1) 


(4) 



(5) 


Equation (2) is valid when the value of the secondary exit stagnation pressure 
is sufficiently high to ensure sonic or supersonic flow at the secondary nozzle 
exit. 
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The normal component of the momentum thrust F is developed from 

n 

the following basic relation 


F 

m 


n 


= p A 
^s s 


' 2 180 ; 


p A cote 

00 S 


( 6 ) 


where e is the injection angle (Fig. 1) measured in degrees. 

Dividing Eq. (6) by the secondary "corrected" mass flow and employing 
basic gasdynamic relationships one obtains the following equation describing 
the normal momentum thrust coefficient 




The forces due to the interaction of the secondary jet with the primary fluid 

are spread in a multiplicity of directions. The resultant force, however, may 

be described as a force acting normal to the primary nozzle wall, F. Work 

in 

performed by Aerotherm personnel (Ref. 19) has shown that the "blast wave" 
theory of Sakauri (Refs. 20 and 21) is valid for a range of "moderately high" 
secondary mass flow rates but that poor agreement with test data was noted 
for the "low" secondary mass flow rates. The work also showed that for 
"low" injection rates the results of Walker and Shandor (Ref. 22) adequately 
predicted the low flow rate regime. 


The interaction thrust coefficient applicable to the high mass flow rate 
situation is derived in Ref. 19 and is herein referred to as Theory 1. 
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The interaction thrust coefficient generated by Theory 1 is given by 



1 




( 8 ) 


where g is an empirically determined correction factor to account for the dif- 
ference between the injection Mach number and the hypersonic Mach number 
for which the blast wave theory was developed. For the freestream Mach 
number range of interest g has been found to be essentially constant, (g =0.7, 
Ref. 23). 

As shown by Dahm (Ref. 24) the value of V is calculated from 


w is calculated from 


where 


V = 1.175 y - 1.05 
P 


, \ 3 / 4 
u = Wj (u 2 ) ' 
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The authors of Ref. 22 have developed a theory (referred to as theory 2) 
to describe the low secondary mass flow rate situation by coupling two-dimen- 
sional linearized supersonic flow theory with the one -dimensional flow theory 
influence coefficients to obtain the resulting expression 
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To determine the range of applicability of Eqs. (8) and (13), Dahm recom- 
mends that one should perform both calculations and conservatively choose the 
smaller value of the interaction force coefficient. 


The only remaining item to consider is the reduction in the calculated 
interaction forces due to the circumferential spreading loss. This empirically 
determined factor is necessary to account for the observed overprediction of side 
forces by both theories 1 and 2. Dahm shows that the factor ( r \ ( j ) ) is reasonably 
well described by 



(15) 


For lack of better data, r) ^ has been assumed to vary with geometry only, 
independent of secondary mass flow rate. 


The previously calculated momentum force coefficients may be resolved 
into components acting parallel (C ) and perpendicular (C ) to the primary 
nozzle centerline, as follows x ^ 



-C m cos(a+e ) + C mn sin(a+e ) 


(16) 


C = C sin(a+e ) + C cos(a+e ) (17) 

rriy m ' mn ' ' 

The forces due to the interaction effect may be characterized by a single 
force acting normal to the primary nozzle wall. Resolving this force into the 
same x and y coordinate system as used in the momentum case, one arrives 
at 


C ix 


C. sina 
in 


( 18 ) 
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C. 

iy 


C in ^ cosa 


(19) 


The final effect of the secondary gas injection may be calculated by com- 
bining the components computed by Eqs. (15) through (18) with the result 


AC 


x 


C + C. 
mx ix 


( 20 ) 


C = C + C. 
y my iy 


( 21 ) 


and the common performance parameter, i.e., the amplification factor (AF) 
obtained from the following relation 


AF 


F /F 
y x 




C m C* 
y s s 

m C* + AC 
P P x 


m 

s 




Figures 2, 3, and 4 show the results of sample calculations for a typical 
secondary injection TVC system. The conditions simulated are: isentropic 
exponent y p = y g = 1.15, secondary injection Mach number = 1.0, primary 
nozzle chamber pressure = 700 psia, primary nozzle area ratio = 7.0, and the 
primary nozzle divergence half-angle = 17.5 degrees. 


The following parameters were varied 


• Secondary mass flow ratio was varied considering a fixed secondary 
injection area with a varying secondary injection pressure. 

• Secondary injection angle was varied from 15 to 75 degrees away 
from the primary nozzle wall. 

• The axial injection location was varied from x/L = 0.25 to x/L 
= 0.75. 
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The calculated data indicate the advantage one gains from injecting 
gas in an upstream direction and also show that moving the injection point 
downstream is advantageous. However, the data at the extreme downstream 
location are questionable due to the incompleteness of the blast wave theory 
(see Ref. 19, p. 2-9). 
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Section 3 

ROCKET NOZZLE FLOW SIMULATION 


3.1 GENERAL 

Future space vehicles presently being investigated will be expected to 
operate from a near vacuum environment through the region of high reentry 
pressures and angles of attack to subsonic flight. To maintain the stability 
and control characteristics of the vehicle over this large range of conditions 
TVC systems and aerodynamic surfaces may be used. 

To design a TVC system capable of meeting these requirements one 
must determine from analytical and empirical means the physical character- 
istics such a TVC system must possess. Section 2 describes some of the 
accepted analytical techniques in use while this section will be devoted to the 
requirements necessary to obtain an adequate TVC test simulation and the 
relationship of these requirements to the proposed MSFC test facility. 
Logically this can be separated into a consideration of the (interrelated) test 
facility geometric aspects and the gas dynamic aspects of the test facility. 

3.2 TVC FACILITY GEOMETRIC REQUIREMENTS 

The primary aerodynamic phenomenon that governs the performance 
of a gaseous injection TVC system is the location and strength of the separated 
flow region generated by injecting a gas into the primary nozzle stream. The 
most important factor pertaining to the strength of the separation region is 
the local Mach number at which the separation occurs. One of the two most 
important parameters affecting the nozzle Mach number distribution is the 
nozzle contour. A typical set of full-scale booster engine contoured nozzle 
coordinates corresponding to a LOX/RP-1 propellant combination is presented 
in Table 5. Using the method of characteristics with real gas chemistry effects 
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accounted for, calculations were performed on this nozzle for a chamber 

5 2 

pressure of 17.237 x 10 N/m (2 50 psia) and an oxidizer-to-fuel ratio of 
2.8 with the resulting internal nozzle Mach number contours being given in 
Fig. 5. Superimposed on these results are the Mach number values obtained 
when method-of-characteristics calculations were made on the proposed 15- 
degree conical nozzle using cold air as the working medium. These results show 
that an adequate simulation of the nozzle internal Mach number distribution in 
the region of TVC injection cannot be obtained by using a conical nozzle to simu- 
late a contoured bell nozzle. Thereafter , it is recommended that a scaled contour 
corresponding to the contour of Table 5 be employed for the MSFC TVC test facility . 

3.3 TVC TEST FACILITY GASDYNAMIC SIMULATION REQUIREMENTS 

In addition to accurately scaling the flight engine geometric character- 
istics, it is necessary to closely duplicate the value of the engine propellant 
ratio of specific heats, y , as this parameter is fundamental to compressible 
flow calculations. Figure 6 shows a plot of the y distribution using LOX/RP-1 
that would occur along the wall of the contoured nozzle of Table 5 and the value 
of the constant y resulting from using cold air as the test medium. Figure 6 
also shows values of y obtained by using a heated variable y simulant gas. In 
this case carbon tetrafluoride (CF^) was selected because it exhibits the de- 
sired low value of y and has good handling characteristics, e.g., low toxicity, 
non-corrosiveness, etc. Figures 5 and 6 show that it is possible to obtain a 
reasonable simulation of the rocket nozzle internal Mach number character- 
istics by duplicating the nozzle contours and selecting a simulant gas whose 
ratio of specific heats is close to the ratio of specific heats of the full-scale 
engine propellants. Figure 7 (Ref. 9) shows the variation of the pu product 
as a function of CF^/N^ mixture ratio for a typical case. 


CF 4 is commercially available from several companies (DuPont, Matheson, 
etc.) and costs about $7.70 per pound. 
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Section 4 

GASEOUS TVC TEST FACILITY REQUIREMENTS 
AND CHARACTERISTICS 


4.1 GENERAL 

This section discusses the mechanical requirements necessary to fabri- 
cate a test facility that will yield pertinent TVC system test data. Suggestions 
regarding the size and materials used for the facility are also given. 

4.2 FACILITY REQUIREMENTS 

To generate meaningful TVC system test data, as previously discussed 
in Section 3.3, it will be necessary to closely approximate the value of the 
ratio of specific heats of the main engine exhaust gas and the TVC system 
injectant gas. Since the propellant combinations for future space vehicles 
are subject to change it is recommended that a means of mixing gases to 
obtain a suitable y be incorporated into the test facility design and that a 
means for heating the primary and injectant gases be provided. By a combi- 
nation of mixed gases heated to the proper temperature, one can very nearly 
duplicate the y value of virtually any propellant. 

Previous test programs investigating gaseous TVC system character- 
istics have shown that the performance of the TVC system is strongly depen- 
dent upon the angle at which the TVC gas is injected into the primary stream. 
Figure 8 (Refs. 26 and 27) presents a typical set of TVC system performance 
data and shows that the level of the side force generated per pound of injectant 
gas is maximum when the injectant gas is inserted at an angle of approxi- 
mately 30 degrees upstream. Therefore, it is recommended that rather than 
inject the gas normal to the nozzle centerline, it be injected in an upstream 
direction at an angle of about 30 degrees. It has also been demonstrated that 
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TVC system performance can be significantly improved by injecting the gas 
through multiple small orifices rather than a few large orifices. Figure 9 
(Refs. 27 and 28 ) shows that increases in the TVC system I from 40 to 80% 
may be anticipated by utilizing multiple orifices. 

When initial consideration was given by MSFC personnel to building a 
TVC system test facility it was tentatively decided to construct it to accom- 
modate a 40:1 scale model of a typical large booster engine and be capable of 

5 Z 

simulating engine thrust chamber pressure levels from 17.237 x 10 N/m 
(250 psia) to 117.21 x 10^ N/m^ (1700 psia). While there is no question that 
such a facility could generate much useful information, further investigation 
disclosed that the capability of the MSFC high pressure air system was unable 
to supply sufficient mass flow to maintain supersonic flow in the nozzle at 
pressures above 41.37 x 10 N/m (600 psia). For this reason and because 
the test facility fabrication and operating costs would be reduced it is recom- 
mended that the model scale be decreased. Figure 10 shows the mass flow rate 
and propellant consumption cost as a function of model size. A throat radius 
of approximately 0.635 cm is suggested as this is still relatively easy to manu- 
facture and would still provide meaningful test data to assess the TVC system 
performance although difficulties may be encountered due to the small in» 
jection port size. Also with regard to model costs it is suggested that con- 
sideration be given to fabricating the nozzle assembly from material that can 
be machined easier than the proposed stainless steel. Aluminum or a filled 
epoxy type of material such as Stycast would be suitable. Also, it is recom- 
mended that the test nozzle be rotatable around its longitudinal axis. This 
would allow the injectant scheme to be modified and would preclude fabricating 
a new nozzle . 

Due to the relatively high temperatures and pressures at which the test- 
ing would be conducted it is recommended that all measurements and all posi- 
tion settings for movable parts such as valves be remotely controllable. Also, 
in order to avert the possible destruction of the test facility in the event of a 
major malfunction it is recommended that a restraining device be provided to 
limit the swiveling of the nozzle. 
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The final observation regarding the TVC test facility is that due to the 
high exit Mach number of the primary nozzle the noise level of the facility 
may be objectionable, necessitating that it be operated only at night. 

4.3 RECOMMENDATIONS 


In summary, the following suggestions and observations are made re- 
garding the MSFC proposed test facility: 


• A means of mixing gases and heating the mixture should be provided 
in both the primary and secondary flow circuits to obtain a suitable 
gamma for testing. 

• The TVC fluid injection angle should be pointed upstream at an angle 
of approximately 30 degrees. 

• A capability for injecting the TVC fluid through multiple small orifices 
rather than one large orifice is desirable. 

• The test facility scale size should be reduced to a throat radius of 
approximately 0.635 cm. 

• Fabricating the nozzle out of aluminum or a filled epoxy material 
should be considered. 

• Making provision for rotating the nozzle on its mounting block will 
give added versatility. 

• Provisions should be made so that all measurements and control can 
be remotely monitored and controlled. 

• A safety ring should be provided to limit total nozzle travel. Testing 
may have to be conducted at night due to possible objectionable noise 
levels. 

Figure 11 schematically shows a facility capable of generating TVC scale 
model test data with the foregoing suggestions incorporated where appropriate. 
It should be noted that the suggested facility is a "blowdown" rig in which pro- 
vision has been made to mix in predetermined proportions a low gamma gas 
from high pressure bottles with MSFC supplied air in an independent plenum- 
heater arrangement. After the gases are bled into the plenums the isolation 
valves shown on Fig. 11 are closed and a suitable amount of time is allowed to 
ensure complete mixing of the gases. After mixing has taken place, the gases 
will be heated to the proper temperature to obtain the desired gamma at which 
time the downstream pressure regulators will be opened and a test run made. 
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Further explanation regarding the facility's operating characteristics 
will now be given. Figure 12 shows two possible ways of preparing the test 
gas to obtain the desired operating conditions (P c = 250 psia, T^. = 600°F). 
Concept 1 consists of a high pressure reservoir containing alumina (Af^O^) 
pebbles. These pebbles have sufficient thermal capacitance and a high enough 
convective heat transfer coefficient to maintain a constant exit temperature of 
the gas. The pebbles are heated in some manner (blowing heated air through 
the system, calrod units, etc. ) to the desired test temperature at which time 
the simulant gas is allowed to flow through the system, down to the pressure 
regulator, and on into the primary portion of the test assembly. Concept 2 
is similar to Concept 1 except that in this case the gas is admitted to the 
reservoir and then heated. See Appendix D for the governing equations and a 
sample calculation for this system. Concept 1 is the recommended design 
primarily because the convective heat transfer coefficient of the pebbles rises 
sharply with the increased flow velocity of system 1 as opposed to system 2 
thereby allowing one to use fewer pebbles and to reduce the total amount of 
energy input to the reservoir. 
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Table 1 

SITVC MATHEMATICAL MODELS 

GAEC (Kaufman, 1968), 2-D 
MDAC (Barnes, et al. 1967) 2-D 
Martin (McDonald & Garbrick) 1966, 2-D 
DVL (Maurer, 1965) 2-D Modified 
U of M (Amick, et al. 1965), 2-D 
GD (Dershin, 1965), 2-D 
ARO (Strike, et al. 1963) 2-D, Circular 
Lockheed (Hair & Baumgartner) 1964, 3-D 
Aerotherm (Dahm, 1967), Circular 
Vidya (Dahm, Mitchell, 1964), Circular 
TRW (Broadwell, 1963) Circular 
NASA (Sterrett, et al. 1966) 2-D 
Lockheed (Carter & Culp, 1970), 3-D 
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Table 2 

BASIC ASSUMPTIONS USED IN LMSC THREE-DIMENSIONAL JI MODEL* 

1. The inviscid flow properties upstream of the separation points (pre- 
stream) are known a priori. 

2. The pre-stream flow is two-dimensional or axisymmetric, including 
the case of local two-dimensionality with respect to a particular 
streamline (pseudo two-dimensional). 

3. The nozzle wall frictional effects on plume properties are presently 
neglected. 

4. The shear layer which characterizes flow separation is of the constant 
pressure type and is either pure laminar (Chapman) or fully developed 
turbulent (Korst). The shear layer developed along the plume boundary 
is considered to be the latter in view of the highly turbulent nature of 
most exhaust plumes. 

5. The average flow velocity within the separation region is very low 
(dead-air) and the average pressure therein can be represented by 
the plateau pressure. Thus, the stream-side and plume-side shear 
layers are subject to the same external pressure (plateau pressure) 
under steady-state conditions. 

6. The Prandtl, Schmidt, and Lewis numbers are unity, thereby render- 
ing the generalized Crocco relation valid. The variation of species 
concentrations or relative mixture ratios across the shear layer can 
be determined in a similar manner. The inclusion of mixing effects 
in the 3-D program has not been completed. 

7. The base temperature and the fuel-to-air ratio can be defined where 
the injector angle is such that the dividing streamlines of the sepa- 
rated flow shear layer and the plume boundary intersection produce 
weak oblique trailing shocks. Although afterburning is not likely to 
occur in the shear layer at high altitudes where long reaction and 
ignition delay times prevail, it is argued that equilibrium chemical 
reactions can be postulated in the dead -air region due to its low flow 
velocity and short characteristic ignition length. 


jJ* 

Lockheed, Carter & Culp, 1970. 
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Table 3 

BASIC INPUT/OUTPUT TO LMSC THREE-DIMENSIONAL SLOT JET PROGRAM 


BASIC INPUT 
Local Mainstream 

e Mach Number 

• Static Pressure 

® Static Temperature 
9 Gas Constant 
e Ratio of Specific Heats 

Exhaust Jet Properties 

• Nozzle Exit Mach Number 

• Stagnation Pressure 

• Total Temperature 
o Gas Constant 

• Ratio of Specific Heats 

Slot Geometry 

® Length 

• Width 

• Injection Angle 

• Area Ratio 

• Nozzle Half Angle 

BASIC OUTPUT 
JIM 

• Upstream JI Force 

• Upstream Amplification Factor 

• Specific Impulse 

• Downstream JI Force 

• Jet Penetration Height 

• Plume Induced Separation Distance 

• Upstream Plateau and Peak Pressures 

• Downstream Plate Pressure 

JIM with Mixing (Current Development) 

• TVC -Propellant Mixture Ratio in Upstream 
Separated 

• Species Concentration 

• Equilibrium Temperature 


34 


LOCKHEED -HUNTSViLLE RESEARCH & ENGINEERING CENTER 






LMSC-HREC D306359 



Local Mainstream Mach Number, M^ 

• A increases with M. 

c 1 

• P^/Pj increases with M^ 

• A^ dependence decreases as M^ increases 

• I increases with M. increasing 

sp 1 

Local Mainstream Static Pressure, P^ 

• A c increases with Pj increasing 

• A c increases with P Q> /P^ decreasing 

• P Q /P decreases upstream separation distance 

‘J . 1 

• P^/P^ decreases 

• Jet penetrates proportionally higher 

Local Mainstream Static Temperature, 

• Independent as long as y^ remains constant 

• Significant with respect to external burning 

Local Mainstream Gas Constant, 

• A^ independent of the mainstream gas constant 

• Significant with respect to external burning 

Local Mainstream Ratio of Specific Heats, y ^ 

• A increases with y, 

c 1 

• A„ increases as Y, increases 

c 1 

Jet Exit Mach Number, Mj 

• A^ varies only 1% as Mj goes from 1 to 2.5 

• At M = 4P < P. JI theory violated 

j oo t 

• F. increases and F. decreases 

J 1 


Most sensitive parameters for two-dimensional Jet Interaction (JI), 
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«T* 

g. Jet Stagnation Pressure, P Q . 

• A c increases as P 0 ^ decreases 

• Smaller penetration — *- reduces separation distance 

h. Jet Total Temperature, T 0 . 

• A c is independent of T 0 - same as Tj (y = const.) 

• I increases as the VTo- 

sp * 

i. Jet Gas Constant, Rj 

• A c is independent of Rj same as R^ 

• I g increases as the ^R ■ (and Vv 

* 

j. Ratio of Specific Heats for the Jet, Y. 

• A decreases as V ■ increases 

c J 

• A„ decreased when Y. increased 

c J 


V fc = a^ ^2/y+l = ^2y/y+lRT 0 .=? V fc decreases as y increases 

m = p. A, V. m decreases as Y. increases 

t t t j 

F = I m F decreases as Y. increases 

s Po J 


F = -m V . + A . ( P. - P_ ) 
J J J °j 


C_ decreases with increase in Y., F = P A. C„ 

F y o. t F 

k* Slot Throat Width, d^ 

• A c increases with a decrease in d fc 

• A decreases 7% with a 100% increase in d 

^ t 

• A decreases 5% with a 100% increase in P 0 

j 

• Increased d and P increases penetration 

t o- 

* o 

i. Jet Injection Angle, cc — 15 to 30 

• A c increases as the jet is transversed to slightly upstream 

• A c decreased as the jet is transversed to slightly downstream 

• 9.j decreases as a - increases forward 

J 


|c 

Most sensitive parameters for two-dimensional Jet Interaction (JI). 
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Table 5 

TYPICAL BOOSTER CONTOURED NOZZLE COORDINATES 
Propellant, LOX/RP-1; P c = 17.237 x 10 5 N/m 2 ; (250 psia); 0/F=2.8 


Local Nozzle Radius, R 
(meters) 

Distance Along Nozzle 
Axis, X (meters) 

0.8636* 

0.0000 

0.8819 

0.1110 

0.8906 

0.1367 

0.9034 

0. 1737 

0.9163 

0.2105 

0.9294 

0.2471 

0.9426 

0.2837 

0.9560 

0.3203 

0.9697 

0.3572 

0.9838 

0.3943 

0.9982 

0.4318 

1.0125 

0.4697 

1.0270 

0.5081 

1.0418 

0.5472 

1.0569 

0.5871 

1.0723 

0.6280 

1.0882 

0.6698 

1.1046 

0.7130 

1.1222 

0.7577 

1.1393 

0.803 5 

1.1570 

0.8510 

1.1752 

0.9002 

1.1940 

0.9514 

1.2134 

1.0047 

1.2338 

1.0605 

1.2554 

1.1192 

1.2769 

1.1802 

1.2990 

1.2441 

1.3220 

1.3114 

1.3457 

1.3822 

1.3717 

1.4577 

1.3974 

1.5367 

1.4237 

1.6202 

1.4510 

1.7087 

1.4800 

1.8033 

1.5101 

1.9041 

1.5403 

2.0111 

1.5715 

2. 1258 

1.6057 

2.2510 

1.7114 

2.6874 

1.7537 

2.8868 

1.8368 

3.3148 

1.8855 

3.5977 

1 . 9222 

3.8420 

1.9430 

3.9877 


•it *** . 

Throat point; tangent point; exit point. 
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Fig. 1 - Component Diagram for Secondary Injection Thrust Vector 
Control Forces (Ref. 19) 
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Secondary Mass Flow Ratio 

TVC Side Force Ratio vs Secondary Mass Flow Ratio 
(X/L = 0.25) 
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Side Flow Ratio (F /F. v ) x 100 
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Secondary Mass Flow Ratio 

Fig. 4 - TVC Side Force Ratio vs Secondary Mass Flow Ratio 
(X/L = 0.75) 
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R/R^, t ~ Nondimensional Radial Coordinate 
' throat 
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Fig. 5 - TVC Nozzle Constant Mach Number Contours 
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Fig. 6 - Variation of the Ratio of Specific Heats for a Typical Booster Engine 
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NOTE: Data shown were generated for 

chamber pressure (P ) of 1000. 
r c 

psia. To use for P c conditions 
other than the above the values 



44 


MLdout fram® 


oj o 



SIDE SPECIFIC IMPULSE, 


SIOC PERFORMANCE AS A FUNCTION OF INJECTION ANCLE 


L.MSG-HREC D306359 





SIDE SPECIFIC IMPULSE 
VS 

THRUST VECTOR ANGLE 


THEORETICAL/EXPERIMENTAL DATA COMPARISON 


PRIMARY PROPELLANTS: LOX/RP-] 



EFFECT OF NUMBER OF ORIFICES ON SIDE SPECIFIC IMPULSE 
DATA FROM CHAPMAN AND GRUNDWALD 



^^8' ^ Effect of Multiple Orifices on TVC System Performance 
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High-Pressure CF^ 
Cylinders 


Fig. 11- Schematic of Thrust Vector Control Facility 
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Concept 1 


Concept 2 


Fig. 12 - TYC Test Facility Heater 
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Appendix A 

MODELING THE THREE-DIMENSIONAL EFFECTS 
OF A SLOT JET FLOW FIELD* 


Lockheed internal research, prepared by R. E. Carter 
and M. F. Culp, Sunnyvale, Calif., 1970. 
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INTRODUCTION 


Preliminary engineering estimates and experiments on flat plates have amply 
demonstrated that the injection of a jet of fluid, either reacting or non -reacting, trans- 
verse to the free stream flow over a surface can produce forces greater than the jet 
thrust alone. Such forces are achieved by separating the boundary layer ahead of the 
jet. The use ot non-reacting transverse jets for force generation has commonly been 
termed jet interaction. Numerous jet interaction tests have been conducted with appli- 
cation to supersonic vehicle reaction control systems and thrust vector control of 
rocket motors. Similar jet interaction phenomena exists when a jet is fired forward 
for thrust termination, when gases escape from a missile stage, and when a highly 
expanded rocket motor plume induces boundary layer separation. 

Jet interaction forces have been mathematically modeled by many investigators but 
their methods are limited to two-dimensional or modified two-dimensional slot flow 
situations. Relatively little attention has been given to the behaviour and structure of 
a three-dimensional model. Unfortunately, the application of a two-dimensional model 
to a three-dimensional reaction control system cannot be made without first determining, 
by testing techniques, several empirical correlation factors. It is therefore believed, 
by the authors of this article, that a basic understanding of the three-dimensional 
effects of the jet interaction phenomenon is required before significant improvements 
to the state-of-the-art technology can be made. This analysis provides a three- 
dimensional model that can be easily expanded to include other jet interaction effects 
such as external burning, viscous mixing, and downstream effects. 
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A* 


“g 

A CD 


a*. 

J 

AR 


b 

b/dt 

BA 


Cn 


C N, A 

C N,R 

dA 

dt 


dy 


dQ s 

FU 

FIFE 

FJ 

FJO 

FN 


SYMBOLS AND NOMENCLATURE 
sonic area 

cross-section area for end flow model 

geometric area at slot span, figure 4 

area where slot span flow reaches freestream conditions 

jet gas velocity (M = 1) 

aspect ratio or b/ dt 

slot half span, figure 2 

aspect ratio - slot half span to width ratio 

width of modified two-dimensional region II, figure 2 

thrust coefficient of jet nozzle 

Maurer amplification factor 

normal force ahead of slot jet 

normal force of jet alone 

surface area over which end flow acts 

width of nozzle throat* figure 2, 3 

increment of span - unit length* figure 4 

slot span Sow 

interaction force in front of the slot 
interaction force at the slot end 
jet reaction thrust of slot 
vacuum jet reaction thrust of slot 
incremental strip force 
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h, or h 


(y> 


“b 

i 

J 

JIB 

K 

S’ 

L 

Mco 

Mi 

m 2 

M e 

“i 

m 

o 

N 

P, 


V P 1 


p . 
oj 


R 2 
**l"2 1 1 

R' 


jet penetration height 

new jet penetration height with slot span flow 

penetration height at slot centerline (ellipse semi -minor axis) 

penetration height at slot span end 

number of selected strips 

constant equal to 0.88, reference 6 

jet interaction with burning 

amplification factor, C °FJO + — * no end etfects 

amplification factor with end effects 

reference length, distance from leading edge of plate to slot 
freestream Mach number, figure 1 
local freestream Mach number, figure 1 

Mach number of inviscid flow above separation region, figure 1 

nozzle exit Mach number, figure 3 

slot span mass flow rate 

slot span mass flow at freestream conditions 

dy - 1, number of strips = 10 

local pressure before separation 

plateau pressure ratio 

pressure behind the jet 

downstream pressure behind the jet 

jet plenum total pressure 

reattachment pressure ratio 

local freestream dynamic pressure 

body radius of slot (cone or cylinder) 

radius as defined in equation 11, figure 2 

6 6 
Reynolds number (laminar < 1 x 10 , turbulent > 5 x 10 ) 
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s 


T 


1 


V 


1 


V 


2 


V 


e 


x.y.z 


X 

X 

e 

a 

6 

a 

°b 

P 

P' 

e 


y> y. 


€ 


£ o 

€ b 
2 -d 


as defined in equation 7 

Local freestream temperature, u lt 

local freestream velocity 

velocity over the penetrating jet 

nozzle velocity at exit to throat sonic velocity 

slot coordinate system, figure 4 

separation distance, figure 2 

new separation distance with slot span flow 

x-axis for blast wave radius, figure 2 

empirical parameter (suggested 1.2) 

cone half angle 

separation ramp angle 

separation flow angle at slot span end 

/ (1 + 7 ~ - m| ) 


empirical parameter (suggested 0.062) 


inclination of nozzle centerline relative to an axis normal 
to the surface, figure 4 


body angle in z,y plane, 


y N-l + ’ 


_L 

2 


N 


- y 


R' 


ratio of specific heats, freestream or jet 

cross flow parameter 

cross flow parameter at slot centerline 

cross flow parameter at end of slot 

two-dimensional flow, figure 2 

local density before separation 

density over the penetration jet 


N-l 


) 


57.3 
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Jet-Freestream Interaction Flow Field 

The phenomenon of a supersonic jet injected into a supersonic freestream has been 
described by many investigators, references 1-9. Figure 1 schematically defines, 
in a profile view, the major features of this jet interaction flow region. The five 
major areas of interest in this figure are: 

• The plume-induced separated flow in the upstream circulatory region (6) that 
increases the reaction force greater than the vacuum thrust of the injectant. 

• The inviscid jet plume flow that is two-dimensional in the case of a slotted 
entrance duct (9), or axi symmetric, in the case of a square port. 

• The penetration mixing region that results from viscous mixing between the 
freestream and the region above the jet normal shock wave (12). 

• The downstream base circulatory region (11) that also effects the force 
amplification of the reaction control system. 

• The downstream recompression region (19) that is dependent upon the 
over-expanded gas (16) in the penetration mixing region. 
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1 



1 . FREESTREAM 

2. LEADING- EDGE SHOCK 

3. INITIAL BOUNDARY LAYER 

4. SEPARATION SHOCKWAVE 

5. SEPARATED BOUNDARY LAYER 

6. UPSTREAM CIRCULATORY REGION 

7. CIRCULATORY STAGNATION POINT 

8. UPSTREAM JET FLOW BOUNDARY 

9. JET ENTRANCE DUCT 

10. DOWNSTREAM JET FLOW BOUNDARY 

11. DOWNSTREAM CIRCULATORY REGION 


12. JET NORMAL SHOCK WAVE 

13. SONIC LINE 

14. MAINSTREAM BOW SHOCK ORIGIN 

15. MAINSTREAM BOW SHOCK 

16. EXPANSION LINES 

17. INTERFERENCE LINES 

18. BOUNDARY STREAMLINE 

19. RECOMPRESSION REGION 

20. TRAILING SHOCK 

21. DOWNSTREAM BOUNDARY LAYER 


o> 


Figure 1 Jet Freestream Interaction 
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Slot Span Flow 

If the jet-freestream interaction is two-dimensional, there will be no flow around the 
end of the jet exit. When flow does occur around the slot ends, the term "slot span 
flow" is used, and a slot aspect ratio effect upon jet interaction performance is 
apparent. Figure 2 schematically illustrates the elliptical line of flow separation 
and four major regions of interest: 

• The two-dimensional region (I), with no slot span flow, can exist near the 
centerline of the slot for high aspect ratio slots with sufficiently low jet 
penetration. However, in some cases this region may not exist at all. 

• The modified two-dimensional region (II), with increasing slot span flow, 
can exist for low aspect ratio slots and sufficiently high jet penetration. 

• The slot span flow is expanded around the end of the slot in region (III) until 
the sonic slot span flow reaches freestream conditions. 

• The downstream region (IV) that realizes viscous mixing between the end 
flow and the downstream region of figure 1. 

Baseline Two-Dimensional Model 

A jet interaction literature search revealed that the two-dimensional analytical method 
developed by Barnes 1 was the best presently available. This mathematical program 
was therefore selected as baseline for improvement since its format could be easily 
modified to account for three-dimensional effects. This mathematical model has the 
ability to calculate the additional force produced by boundary layer separation that 
results from a gas jet injected through an infinitely long slot in a fiat plate surface 
into a uniform supersonic stream. 
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LEGEND 

T 2-d REGION - MDAC 
AND OTHERS 

n MODIFIED 2-d REGION 

m END EFFECT 

TTT DOWNSTREAM REGION- 
HIGHLY EMPIRICAL 
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Figure 2 View Looking Down on Jet Transverse to Flow 
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The baseline model uses a two-dimensional momentum balance between a transverse 
slot jet and the surrounding supersonic freestream. The momentum balance consists 
of equating the drag on the equivalent body formed by the transverse jet interaction to 
the change in axial momentum of the transverse jet. This model assumes that the 
viscous portion of the drag is proportional to the reattachment pressure coefficient 
while the inviscid contribution (deflection of the freestream) is proportional to the 
plateau pressure coefficient. The total drag coefficient is related to these parameters 
by empirical constants. The model computes the jet penetration height (h Q ) of the 
injected gas, the point of freestream flow separation (X g ) in front of a flat plate slot 
with end plates, and the force amplification (K). These parameters are a function of 
the plume-induced separated flow for local laminar or turbulent freestream Reynolds 
numbers, jet and freestream gas properties, and freestream Mach number. 

Jet Penetration Height 

The first modification to the two-dimensional baseline model involved a new 
method 2 of calculating the penetration height (h) of a transverse jet. To do this, the 
plateau pressure (P ) can be computed using the curve fit of Barnes 1 for pressure 
coefficient. Both the laminar and turbulent flow values have been mechanized and 
can be computed as a function of Mach number. The plateau pressure is modified 
by a new constant (£) which is representative of the reattachment pressure on a 
forward facing step. The pressures on the back side of the jet are analogous to the 
base pressure of a step and the constant (a) is used to represent the base pressure to 
freestream pressure ratio. The sum of the forward and backward face pressures is 
then the total force acting on the jet in the streamwise direction. 
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The penetration height (h) calculation is therefore no longer dependent on the baseline 
equivalent body drag coefficients and empirical constants of Barnes 1 : 



2 

It is suggested that values of ft • = 0. 062 and c* = 1. 2 give good agreement with jet 
interaction force data. A study was conducted to determine the sensitivity of the jet 
penetration height (and hence the boundary layer separation distance) to the emperical 
parameter, /? The sensitivity of the computed boundary layer separation distance 
to different values of /? ' for different nozzle configurations, and test data 3 , is pre- 
sented in figure 3. 

It was also found that the value of /?' = 0.062 gives adequate agreement with centerline 
reattachment pressure of Maurer 4 . The reattachment pressure 2 is: 
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JET TOTAL PftESSURE/AMBIENT PRESSURE, P ./ P, 
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Figure 3 Boundary Layer Separation Distance vs. Pressure Ratio 
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Three-Dimensional Flow Field Model 


When flow begins around the slot ends, slot span flow occurs, and an aspect ratio ef- 
fect upon jet interaction performance is apparent. Treatment of slot span flow was 

4 

first suggested by Maurer by integrating the flow from the centerline of the slot to 
the slot end, as shown in figure 4. To start the computation procedure, a separation 
distance (X /L) at the centerline (distance the flow separates ahead of the jet) for two- 
dimensional flow is computed. This distance depends on the local flow, laminar or 
turbulent, Reynolds number (R^), the Mach number (M^), the slot half span to width 
ratio (b/dt), the ratio of specific heat ( y ), and the given jet pressure ratio P^/P^). 
Estimates of slot span flow are made using the techniques of Maurer, and new separa- 
tion distances (X) are computed along the slot span. Maurer's oil-film photographs 4 
show, except for the case of a large aspect ratio slot and small jet pressure ratios, 
the separation line to be elliptic. In this case, the nozzle slot itself is parallel to and 
in the vicinity of the major axis of the ellipse. The jet penetration height (h ), figure 

y 

4, is therefore assumed to be elliptical when slot span flow is present. 


A cross flow parameter. 


€ = 


dQ s 

h (y) P l V l * 


OS £s 1.0 


(5) 


specifies a fraction of mass flow dy h^ ^ V 1 entering the control volume surface (1), 

figure 4, which does not exit through surface (2) after deflection at the separated flow 

wedge, but spills out the sides of the control volume as dQ , figure 4. Thus, in each 

s 
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section dy, there exists a cross-stream flux fed by a part of the main stream which 
entered the control volupie. 

All freestream flow along the slot centerline shown in figure 4 is initially assumed to be 
two-dimensional (cross flow = 0). For the three-dimensional solution, the cross 
flow parameter ^ b is determined at the end of the slot span such that the penetration 
ellipse equation is satisfied. An iterative solution is used for the ellipse equation 


1 - ( hj/1^) 2 - 1/(1 + h b /(b tana)) 2 


( 6 ) 


such that for the given value of € q the penetration height of the slot centerline is cal- 
culated. A search is then conducted to determine e b which solves the above equation. 
For each guessed value of e b the penetration height at the slot span (h b ) and the sep- 
arated flow or wedge angle ( <7 b ) is calculated. The program then increments the cross 

flow parameter ( e )between the initial value at the slot centerline ( f =0) and the 

o ' 

value determined at the slot end ( e b ). The details of this step increment are such 
that each value of € is related to a distance along the slot span (y) by the penetration 
ellipse equation. Thus, relatively large dy increments will exist near the center of 
the flow field where the slot span flow is negligible (approximately two-dimensional) 
and small where the cross flow is significant. 
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The difference between the freestream mass flux entering the control volume surface 
(1), figure 4, and the mass flux over the penetrating jet through surface (2) must 
equal the lateral flow component dQ g . 

Thus - «*>1 V l h o ' fl. v 2 W COS ° P) 

Tbe ratio i> s /h o may be derived from equation (7) and a momentum balance across 

the oblique separation shock, thus 


...frg - 

h 


1 



( 8 ) 


Finally, the product of equations (1) and (8) provides the new penetration height (hs) 
that when divided by Tan a will define a new separation distance in front of the slot (X ) 
The two-dimensional separation and penetration distance is then modified by the lateral 
cross flow at N along the slot span. In all cases, the slot span flow is redefined until 
a sonic slot span flow condition exists at the end of the slot. This is accomplished by 
the computer in the following manner. The geometric area (Ag), figure 4, formed by 
the penetration height (h^) and the separated flow distance in front of the slot (X^) is 
checked against the choked flow area (A*) associated with the fraction of mass diverted 
back into the separated region. 
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A - h b X b <9) 

8 2 


If A equals A* at the slot span, then the three-dimensional solution has been obtained. 
6 

If A* is less than A , then the cross flow parameter at the centerline (e Q ) is increased 

s 

until A* equals A . If A* is greater than A , a two-dimensional flow exists in the 
g g 

center section of the slot span and a new slot semi -span (BA) (see figure 2)is computed 
whose centerline value is = 0. 

End Flow Models 


The proposed mathematical model treated the flow ahead of the slot from its center 
line out to the end of the slot. Sections I and II, (see figure 2). That portion of the 
freestream that does not pass over the slot jet is diverted outbound in front of the 
slot, creating slot span flow. Region III (see figure 2) remains to be defined and 
extends beyond the end of the slot. 
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Maurer^ proposes swinging a radius equal to the separation distance at the slot span 
for estimates of forces beyond the slot end. The end flow equation in the computer 
parametric studies used this technique and is as follows: 


FUE = 


7TL 

4b 


/M 2 / p 2 - P 1 \ 
\ L / \ / 


( 10 ) 


14 

Little is known about the end flow of a slot * and it is not intended at this time to 
precisely describe the flow, but to propose a different model suitable for refinement 
when experimental data becomes available. Nunn and Karamcheti propose a blast 
wave solution for a round jet on a flat plate or on the side of a nozzle. In this analysis 
the slot span flow at the end of the slot is analogous to the jet on a plate or the side of 
a rocket nozzle, and the blast wave method defines the slot span flow penetration into 
the freestream. This penetration distance (R e ) when traced on the surface adjacent, 
to the slot end, is then the separation distance with respect to the slot end. 


R = 
e 




y m. X 1 /* 

J. L Aocl 

J m / 

o / 


hr 


<ii) 


The radius is then computed along the side of the slot to the point where the slot span 
flow (dQ ) has been expanded insentropically from sonic conditions at the end of the slot 
to the freestream Mach number (see figure 4). The pressure in Region HI, 

Figure 4, is assumed to be the plateau pressure (P^) Region n at the slot end, 
decaying to freestream at the line of boundary layer separation. The pressure 
distribution and resulting cross-sectional area (A ) is assumed to be that of the 


17 


LOCKHEED MISSILES & SPACE COMPANY 



LMSC-HREC D3C63 59 


extended penetration ellipse. The following integration is then made over the surface area 
beyond the slot end for the force: 


FUE 



< P 2 ~ P l> dA 
q x bL 


( 12 ) 


which, when pressure-area terms are included and integration completed, is: 


FUE 



L 

tan o 


1- <V X e> 5/2 


cos 6 cos 0 


(13) 


The above technique will allow adjustment of the pressures over the areas affected, 
the shape of the penetration ellipse, and provide a first real approximation of the 
flow separation distance beyond the end of the slot. 

Surface Curvature Approximation 

The necessity of considering the effects of surface curvature on amplification factors 
are the result of practical applications of a jet steering system. Nunn 5 estimates 
the effects of body curvature for orifice amplification factors using a blast wave 
analogy to the problem. A method has been devised here for the slot jet with slot 
span flow. The analysis is proposed as a method for extending the effects of body 
curvature from an orifice to slot jets with slot span flow. 
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The total force (FU) upstream of the slot is derived by integrating the flat plate 
component force (FN) for each strip shown in figure 4. For conical or cylindrical 
bodies, the average component force normal to the x, y plane is computed prior to 
integration as follows: 


i 

FU ^ FN cos d Cos 6 
N=1 


(14) 


The force (FUE) at the end of the slot may be corrected in a similar manner for body 

5 

curvature by using the technique proposed by Nunn. 

Finally, the total amplification factor (K T ) is computed and represents a vector 
normal to the x, y plane. 


_ FU + FUE + FJ Cos 0 Cos d Cos 9 
T FJO 


(15) 
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Comparison to Maurer's Test Data 
4 

The experimental data of Maurer has been replotted and the amplification factors 
adjusted to use slot jet coefficients referenced to a vacuum. Maurer is not exactly 
clear about his reference thrust. On the one hand, he refers to the reference slot 
Jet as one coming from a real nozzle, but also implies that his amplification factors 
(C Q ) are referenced to nozzle characteristics of a jet exhausting into an atmosphere 
rather than a vacuum. The following equation was derived from ideal nozzle equations 
and used to adjust Maurer's experimental amplification factors to a vacuum reference. 


C (Maurer) y + 1/y 

n * = _1 L (16) 

K (Vac Ref) V*. 


The nozzle exit velocity to throat velocity, V /a* , is a function of jet chamber 

e J 

7 

pressure and can be read directly from compressible flow tables , with this, the 
conversion is easily accomplished. 


Figure 5 shows the comparison of Maurer's test data (symbols) and the computer 
program results at the centerline with no end effect. Maurer's test data are based 
on centerline integrated pressures for amplification factors, and because of this, no 
end effects can be evaulated using his experimental data. The tendency to under predict 
the amplification factors at high jet pressure ratios is evident. Further refinement of 
the mathematical model to better predict separation pressures at low Mach numbers 
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Figure 5 Maurer, Mach 1. 57 Comparisons 
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should improve this trend. When applying a jet interaction control system to a real 
time trajectory, it was found that low Mach numbers are usually associated with low 
altitude flight and low chamber pressure ratios. In light of this fact, the dropoff of 
predicted amplification factors at high jet pressure ratios may not be as important to 
control system studies as it may first appear. 


Comparison To Romeo's Test Data 

g 

Romeo's experimental data were ootained for a slot jet transverse to a Mach 6 
freestream flow on a flat plate. Pressure taps were located ahead, to the side of, and 
in back of the slot. These pressures were integrated over the plate and compared 
to the thrust of the jet with no flow. Amplification factors are not given, however, 
the normal force coefficient for the integrated pressures ahead of the jet and the 
normal force coefficient of the jet were shown. The following equation was derived 
from the ideal gas equations for a sonic nozzle to convert the data of Romeo to ampli- 
fication factors: 



Figure 6 presents a comparison of the present analysis with the test data of 

Romeo. A reasonable correlation was found for the variables investigated. A closer 

correlation may have been found by adjusting the high Mach number constants ( 0 or a ) 

used in the separation characteristics. The sensitivity of one of these factors was 

2 

explored and presented in the Spaid and Zukoski model for separation. 
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Figure 6 Romeo, Mach 6 Comparisons 
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Speculation of test equipment problems in the experiments that Romeo conducted 
comes from the dramatic rise and fall in the amplification factor for small throat 
widths (dt). Sterrett** also noted this effect. It is probable that for the small slot 
widths and high jet pressures, the slot was significantly distorted. If this occurred, 
the jet may have been tilted forward into the stream with an unknown increased slot 
width and basic jet thrust change. The construction of the slot in the test rig was such 
that the back edge of the slot was the edge of a thin flat plate. When this plate bulges 
due to high jet chamber pressures, the slot would appear to be aimed forward (see 
Reference 8 for test rig details). Then, as demonstrated by Maurer 4 , where the 
slot was aimed into the stream, the amplification factor increased substantially. 
Because of this uncertainty, the test data at throat widths of 5x10 inches or less 
are not shown. 


Jet Interaction Computer Studies 

The jet interaction computer program used for this study uses the method of Barnes 
combined with Maurer's slot span flow and Spaid and Zukoski separation techniques 
as described in the preceding comparisons. End effects were estimated using a 
radius at the slot end. No computer capacity limitations were experienced during any 
of the studies, including the parametric studies, although many thousands of data 
points were retained for cross plotting. 

The computer run parametric study has been included for demonstration purposes. It 
is typical of a slot jet application to a low altitude (atmospheric) supersonic missile. 
Two types of information were plotted, (1) data used for engineering evaluation of 
the method employed, and (2) data for systems evaluation when conducting preliminary 
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design studies. The following were the basic input values for this parametric 
study: 


Slot jet assumed to be installed on a flat surface. 
y 1.4 freestream 

y j 1.25 nozzle 


* 

a. 


1. 805 velocity at exit of nozzle to throat 



_b 

dt 

dt 

L 


b_ 

L 


0 

R = lxlO 7 
n 


40-100 chamber pressure to freestream static 
2 to 6. 6 

10 to 90 aspect ratio 

0. 00288 throat diameter to reference length 

0.00144—0. 131 half span to reference length 
900 °R freestream temperature 
39° nozzle pointed forward from vertical 
Local Reynolds number at nozzle, turbulent 
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Figures 7 thru 10 present the results of the parametric study for engineering use. 

Amplification factors for the slot span are integrated values along the span with a 

quarter circle radius at the ends for end effects . Centerline amplification factors 

are shown by assuming centerline conditions existing uniformly over the slot span 

with a radius at the slot end for end effects. For information, two-dimensional 

values with no end effects are shown. X g is the separation distance ahead of the 

slot referenced to the distance from the origin of the boundary layer to the slot. 

h A is the ratio of the separation height at the slot (with slot span flow) to the 

separation height if no slot span flow existed. P /p is the ratio of reattachment 

K 1 

pressure to local freestream pressure ratio. 


Figure 11 is a plot of amplification factors for a slot with end effects as a function 
of aspect ratio, Mach number, and chamber pressure. The forces a slot jet system 
can produce with respect to a slot exhausting to a vacuum can be computed using the 
amplification factors presented. For example, 7 = 1.4, P^A^ = 80, Mach num- 
ber = 4 and b/dt of 10, Figure 11 shows K^, equal to 2.25. The control force is 

then FU + FJO = 2.25 C, . P . . dt, where: 

fo oj 




1.268 


(18) 


From the figures, increased aspect ratio and Mach number show increases in 
amplification factor. Conversely, as jet pressure ratios increased, amplification 
factors decreased. The time required on the IBM 1108 computer to conduct this 
study was five minutes including the plotting. A new parametric can be set up for 
the computer in about two hours . 


26 


LOCKHEED MISSILES & SPACE COMPANY 



LMSC-HREC D306359 


— I LEGEND — 

— — SLOT SPAN 
— — — CENTER LINE 
X— X— -X K-2D ~ 

l I 

MACH NO .= 5 
b/L = .0006 

ASPECT RATIO b/d + = 10 






LMSC-HREC D3063 59 



Figure 9 Separation Distance Vs Jet Pressure Ratio 



Figure 10 Reattachment Pressure Vs Jet Pressure Ratio 
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Conclusions 

• A large number of jet interaction variables were evaluated and the amplification 
factors computed show reasonable agreement with experimental data. 

• The empirical constants, ft = 0.062 and a = 1.2, of equations (2), (3) 
were found to provide adequate agreement with centerline pressures and 
boundary layer separation distance estimates. These constants were then 
used in the basic momentum balance to determine the jet penetration height. 

• Slot span flow greatly affects the prediction of amplification factor for low 
aspect ratio slots at a significantly high jet-to-freestream pressure ratio. 

• Slot span end flow greatly increases the prediction of total amplification factor. 
The end flow surface area has been approximated by swin g in g a radius (quarter 
circle) equal to the separation distance at the end of the slot. The blast wave 
theory is proposed for the small aspect ratio slots at significantly high pressure 
ratios. 

• Total amplification factor is reduced by body curvature. Predicted amplifica- 
tion factors can be estimated by correcting the differential force increment 
for the warped surface prior to integration. 

• The three-dimensional model developed enables future parametric study to 

include the effects of viscous mixing in the separation and penetration miying 

regions, and provides the flow conditions necessary for the prediction of down- 
stream effects. 
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Recommendations 

• Incorporate a viscous modeling technique such that the amount of freestream-jet 
miving in the separated region can be estimated. 

• Model the downstream flow field specifically to determine the reattachment region 
distance from the slot. 

• Complete a study of pressure distribution on the immediate aft side of the jet 
using a method dependent on freestream and flow conditions above the jet. 

• Define a technique and estimate the effects of jet interaction with burning (JIB) 
on the amplification factor. Both the upstream and downstream regions should 
be considered using a reactive gas jet. 
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Appendix B 


A detailed description of the program input instructions follows: 

• Detailed Input Guide 

• Sample Program Input 

• Sample Program Output. 


Card 1 


Format 7E 10. 5 


Column 

Parameter 

Definition 

1-10 

SPR 

secondary-to-primary total 
pressure ratio 

11-20 

XL 

ratio of distance from primary 
nozzle throat to the injection 
point divided by nozzle length 

21-30 

EPS 

secondary injection angle (deg) 

31-40 

SM 

secondary injection Mach number 

41-50 

CMS 

secondary nozzle flow coefficient 

51-60 

VNS 

secondary vacuum momentum 
specific impulse efficiency 

primary total pressure (psfa) 

61-70 

p T p 


Card 2 

Format 7E 10.5 


1-10 

ARATIO 

primary nozzle area ratio 

11-20 

GAMS 

ratio of specific heats for the 
secondary gas 

21-30 

VNP 

primary vacuum momentum 
specific impulse efficiency 
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Column 

31-40 

41-50 

51-60 

Card 3 

Format 

1-10 

11-20 

21-30 

41-50 

Card 4 

Format 

1-10 

11-20 

21-30 

31-40 

41-50 

Card 5 

Format 

1-10 

11-20 

21-30 

31-40 


Paramete r 

CDP 

SMI 

PSL 


Description 

primary nozzle flow coefficient 
initial secondary mass flow ratio 
local static pressure (psfa) 


6E10.5 


WMP 

primary gas molecular weight 

WMS 

secondary gas molecular weight 

TTP 

primary gas total temperature ( R) 

WPRI 

primary mass flow rate 


7E10.5 


TTS 

GAMP 

PM 

S 

AS 


secondary gas molecular weight (°R) 

ratio of the specific heats of the pri- 
mary gas 

local primary flow Mach number 

interaction distance (ft) 

/ 2 . 

secondary flow area (ft ) 


5E10.5 


CFP 

ALPHA 

SMINC 

SMMAX 


primary nozzle thrust coefficient 
nozzle divergence angle (deg) 
secondary mass flow increment 
maximum secondary mass flow 
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SAMPLE PROGRAM INPUT 


■ill) A I A . 1 L 

C HU L fINI 







OAlA HROCESS^ l t . EL 

3 






OOUU J 1 

iqlM». H , 1 L A 

.A 






U0lHJLi2 

ue s>rs»*MSFci. 






U GuL U 3 

itfXUT A 







U G U U C *1 

.V 

.25 

75. 

1. 

; .’3 

.965 / 

100800. 

OGUUOS 

7. 

1.15 

.9768 

.99 

• 02 

13225. 


00UGU6 . 

20. 

20. 

1500. 

100. 




0 U u U u 7 

1 bOO . 

1.15 

2.011 

.326 

•000729626 



OOUUOB 

1 .St 

15.0 

.02 

.29 




UUUULV 

iai X W T A 







Otluu 1 U 

.9 

.25 

60. 

1 9 

.93 

• 965 

100800. 

00U0 1 1 

7. 

1.15 

.9768 

.99 

• 02 

13225. 


OOUU 1 2 

20. 

20. 

1500. 

100. 




UOUU 1 3 

I SCO. 

1.15 

2.011 

.326 

•000729626 



OOUG 1 H 

1 .SI 

: 5.c 

.02 

.29 




UUUU 1 b 

ml* UT A 







OOUCi 1 o 

.V 

• 25 

95. 

1 • 

.93 

.965 

100800. 

GGUU 1 7 

7. 

1.15 

• 9768 

.99 

.02 

13225. 


uouu 1 a 

20. 

20. 

1S00. 

100. 




oouo i * 

1 bOO . 

1.15 

2.011 

.326 

•000729626 



UOu02U 

1 .SI 

15.0 

.02 

.29 




U U u b i 1 

w X VI T A 







J0U022 

.9 

.5 

7S. 

1 . 

.93 

.965 

100800. 

UL) U U 2 3 

7. 

1.15 

.9768 

.99 

.02 

6350. 


U0UU21 

20. 

20. 

1500. 

100. 




UUUU2 b 

1 SOU. 

1.15 

2.905 

.2171 

.000729626 



0 OU U 2 o 

1 .51 

15.0 

.02 

.29 




OOuU 2 7 

WXWT A 







UUU U 2 0 

.9 

.5 

60. 

1 . 

.93 

.965 

100800. 

UUuOi’V 

7. 

1.15 

.9768 

.99 

.02 

6350. 

UCGU3U 

20. 

20. 

1500. 

100. 




J U c 0 J 1 

1 5 30. 

1 . is 

2.905 

• 2171 

.000729626 



UOULi 3 2 

1 .51 

15. U 

.02 

.29 




U U U U 3 3 

I3XLT A 







UuuiJ J *i 

.9 

• 5 

95. 

1 . 

.93 

.965 

100800. 

UOuG 3b 

7. 

1.15 

• 9768 

.99 

• 02 

6350. 


Li 0 b U 3 6 

-0. 

20. 

1500. 

100. 




GUU03 > 

1 500. 

1 . IS 

2.905 

• 2171 

• 00u7 29626 



J 0 U Li 3 0 

1 .51 

15.0 

.02 

.29 




LiOuU 3 V 

■3X L T A 







UUuUMU 

.9 

.75 

75. 

1 . 

.93 

• 965 

100800. 

OUuuH 1 

7, 

1.15 

.9768 

.99 

• 02 

365.6 


l, o o a *i 2 

20. 

20. 

1500. 

100. 




(J 0 U U 4 3 

1 SOU. 

1 . 15 

2.680 

. 1083 

•000729626 



atijuaa 

1.5 1 

15.0 

.02 

.29 




UO uO M *> 

v X w T A 







uuul)‘lt> 

. 9 

.75 

60. 

1 1 

.93 

• 965 

100800. 

OUuU ‘i / 

7. 

1.15 

.9768 

.99 

• 02 

365.6 


UU uU *« d 

2u • 

20. 

1500. 

100. 




JOUU'iV 

1 5CU. 

1.15 

2.680 

.1063 

•000729626 



U'UUOijG 

1 .51 

15.0 

.02 

.29 

... 

*' 


J G L G g l 

•» X T A 







JOuG j2 

. 9 

• 7S 

95, 

1. 

.93 

.965 

100800. 

iJQuU'i 3 

7. 

i . : s 

.9768 

.99 

• 02 

365.6 


i j 0 u 0 H 

20. - 

20. 

1S00. 

too. 




UOuUbb 
UOuU b 6 

0. 

i • b 1 

1.15 

15.0 

2 .680 
.02 

.1083 

.29 

•000729626 

— 
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. AS 

CFMP 

CFMN 

thet 

smr 

7 • 2 96 26-0*1 

1 « 004 4 3 + 00 

6,75623*02 

1.39806*00 

2,00000-02 

cfini 

CFIN2 

CF1N 

AF 

PTS 

2.84115*01 

8.75706"01 

8,75706-01 

\ 

I , 22026*00 

4,67742*04 


cfy 

C F MX 

cfmy 

ALPHA i 

.... 

3. 14210-01 

1 .85029 + 00 

6.75612*02 

1 ,QQH43t00 

1,50000*01 


FS/FA 

2.44056-02 . 


w 

I 

4^ 


AS 

CFMP 

CFMN 

thet 

smr 

7 • 29626-0** 

1 .14275*00 

1,24627-01 

2,96773*00 

4.00000-02 

CFINI 

CF1N2 

CFJN 

AF 

PTS 

1 ,42057*01 

8, 75706-01 

8,75706*01 

1,30483*00 ’ 

9,75464*04 

CFX 

CFY 

CFMX 

CFMY 

ALPHA 

3.51275-01 

1,98862*00 

1,24626-01 

1.14275*00 

1 .50000*01 


FS/FA 

S. 21930-02 


W 

> 


u 

r 

w 

*d 

JO 

o 

§ 

> 


S 

o 

c 


H 

c! 


H 


r 

2 

c/i 


ffi 

jo 

M 

O 

d 

OJ 

o 

O' 

U) 

cn 
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*Fon,lS HA 1 1«. N 

'i 5 l uoY-12/1-- . ! u 6 ( , 0 ) 


MAIN PROGRAM 

STORAGE U>E U I CODE. I I I 0012511 DATA(O) 0002421 BLANK cOMMQN(2) 000000 


0 

1 

►— 

1 


EXTERNAL REFERENCES (BLOCK," NAME) 


0u03 

NINTHS 

CL OR 

NRDOS 

LOOS 

N I 0 2 S 

0006 . 

NO X P 6 S 

LLD7 

S * ft T 

C o 1 0 

T a n 

Do 1 1 

COS 

00 1 2 

SIN 

0013 

A Tan 

00 1 4 

N..DUS 

00 1 5 

NSTOPS 


storage assignment (Block, ttpe, relatjve location, name) 


0000 


000 1 62 

10F 

0001 


00007 1 

IDOL 

0001 


001245 

1 0 1 L 

0001 


000325 

1 1 IL 

0001 


000400 

1 12L 

uuOO 


OUOl 64 

12F 

0000 


00016* 

IMF 

0000 


000170 

16F 

OOOO 


000172 

10F 

oooo 


0001 16 

200F 

0000 


OoOl 30 

20 1 F 

0000 


000133 

202F 

0000 


000145 

2q3F 

oooo 


000157 

204F 

oooo 

R 

000060 

A 

uuOO 

K 

0001 1 2 

AF 

0000 

K 

000030 

alpha 

0000 

R 

000053 

AP 

oooo 

R 

000010 

ARATIO 

oooo 

R 

000026 

AS 

UoUO 

R 

000047 

ASEFF 

0000 

H 

OOOuSR 

ASTARP 

0000 

R 

000U6 1 

B 

oooo 

R 

000062 

C 

oooo 

R 

000013 

CDp 

0000 

R 

QUO 102 

CFIN 

0000 

R 

OOUU7S 

CF1N1 

0000 

R 

OOOlUO 

CFIN2 

oooo 

R 

000103 

CFIX 

oooo 

R 

000104 

CF 1 Y 

0 J .00 

R 

OLOUA 3 

cfmn 

0000 

R 

0000S7 

CFMP 

OOoO 

R 

000 1 06 

CFMX 

oooo 

R 

000107 

CFMY 

oooo 

R 

000027 

CFP 

UOUO 

R 

occur 1 

CFVO 

ouoo 

R 

0001 10 

CF A 

0000 

R 

0001 1 t 

CF Y 

oooo 

R 

000036 

CGAMP 

oooo 

R 

000035 

CGAMS 

0000 

R 

OOLOOS 

CMS 

0000 

R 

000037 

CSTAKP 

OOOO 

R 

□ 0 0 u 4 0 

CSTARs 

oooo 

R 

000056 

D 

oooo 

R 

000055 

OP 

0000 

H 

0 u 0 0 6 s 

E 

UOOu 

R 

000003 

EPS 

OOOO 

R 

000066 

F 

oooo 

R 

0000*7 

<a 

oooo 

R 

000023 

GAMP 

0000 

k 

00001 l 

Gams 

0000 

H 

0001 1*» 

G I MANG 

oooo 

R 

00007 I 

H 

ouoo 

I 

000034 

ICOUNT 

oooo 

R 

oooono 

L 

OOoO 

K 

OOOC 7 2 

p 

0000 

R 

000101 

PED 

oouo 

R 

000024 

PM 

oooo 

R 

000015 

PSL 

oooo 

R 

000007 

PTP 

0000 

R 

QU00H3 

PTS ...... 

0000 

R 

C0007 6 

9 ....... 

oooo 

R 

000105 

R 

oouo 

R 

000025 

S 

oooo 

R 

000004 

SM 

ULUO 

k 

OloO 1 4 

SM) 

0000 

R 

00003 1 

SmINC 

oooo 

R 

000032 

SMMAX 

ouoo 

R 

000033 

SMR 

oooo 

R 

noooo 1 

SPR 

LOUO 

R 

000 1 1 3 

T 

0000 

R 

000044 

TERMl 

oouo 

R 

00 004 5 

TERM2 

oooo 

R 

000046 

terms 

oooo 

R 

000050 

terms 

CUuO 

K 

OOOCS l 

T ERMS 

0000 

R 

0 00 u S 2 

TERM6 

ouoo 

R 

0001 15 

thetag 

oooo 

R 

000020 

T TP 

oooo 

R 

000022 

TTS 

0 JOO 

h 

OoOOAR 

TU 

0000 

R 

00001 Z 

V NP 
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R 

000006 

VNS 
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R 

000074 

A 
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R 

OOOO 1 6 
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0000 

R 

UuOO 1 7 

WMS 

0000 

K 

000021 

APR 1 

oooo 

R 

000042 
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□ 000 

R 

000070 

At 
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R 

000073 

*2 

uO 

ft 

UuGu7 7 

*3 

0000 

R 

uooooz 

XL 














i* c 

2 • C 

3* C 

s» c 
s» c 
*> • c 

7 • C 


• •• INPUT PARAMETERS , 

af ■ amplification factor 

ALPHA ■ PRIMARY NOZZLE LOCAL hall ANGLE (DEGREES) 
AH AT 1 0» PRIMARY NOZZLE AREA RATIO 

AS » SECONOAKY GEOMETRICAL FLO* AREA (S<JUa«E FEET) 

aseff»effect I ve secondary floaahea iswoare feeti 
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JOIUu 

8* 


001 OU 

9" 


0010U 

- .. 10* 

C 

0U10J 

1 1" 

c 

OUIUO 

12" 

c 

UUlOo 

13* 

c 

UU 1 OU 

|9. 

c 

001 00 

1 S* 

c 

OOIOU 

16" 

c 

UO 1 Ou 

17" 

c 

UU 1 OU 

18" 

c 

UO 1 UU 

19" 

c 

UU1GU 

20" 

c 

UO l UO 

21" 

c 

UUIOO 

22" 

c 

UO 1 UU 

23* 

c 

UO 100 

29* 

c 

UU 100 

2b* 

c 

uOIOO 

26* 

c 

OUlGu 

27* 

c 

UU 1 Ou 

28. 

c 

UUIOO 

29* 

c 

00 1 Gu 

30* 

c 

uO 1 Ou 

31* 

c 

uO l 00 

32* 

c 

OOIOU 

3 3* 

t 

uO 1 UU 

39* 

c 

OOIUO 

3 S • 

c 

UO 1 GU 

36* 

c 

uO 1 UO 

37. 

c 

UUIOO 

30* 

c 

UU I UU 

39. 

c 

UO 1 Ou 

9 0 * 

c 

uO I OU 

9 l • 

c 

UO I Uu 

92* 

c 

UU l Ub 

93* 

c 

OOIUO 

99* 

c 

UO 1 Uu 

9 S • 

c 

UU l uu 

96* 

c 

00 i Ou 

97* 

c 

UU I.Ou 

98* 

c 

UO 1 uu 

99* 

c 

00 1 Ou 

SO* 

c 

UO 1 U 0 

SI • 

c 

u o i o 

b 2 • 

c 

00 1 U 

b 3 * 

c 

jo 1 . 

SR* 

c 

j 0 i - . 

S S • 

c 

U U I U 1.. 

b 6 * 

c 

UU 1 ou 

b 7 » 

c 

UO 1 uu 

b fa * 

c 

JOT 0 1 

bV* 


UUIOO 

60* 


JOi 1 R 

6 l » 


uUIZR 

62* 


JUl 32 

63* 



astarPbpnimary nozzle throat area (Square feet> 

COP • PRIMARY NOZZLE FLO* COEFFICIENT 

CFINb MINIMUM VALUE OF CFINl OR CF I N2 _ 

cfini ■ huh secondary massflow induced force coefficient 

CFJN2 ■ LOtt SECONDARY MASSFLOW INDUCED FORCE COEFFICIENT 
CFMN. NORMAL SECONDARY INJECTION MOMENTUM THRUST COEFFICIENT 
CFMPa PARALLEL SECONDARY INJECTION MOMENTUM THRUST COEFFICIENT 
CFMX "MOMENTUM THRUST COEF. PARALLEL TO PRIMARY NOZZLE AXIS 

CFMY "MOMENTUM THRUST COEF. NORMAL TO P R I M *R Y NOZ Z LE AXIS 

CFP - PRIMARY THRUST COEFFICIENT (nO TVC EFFECT) 

CFVO. isenthopic thrust coefficient 

cfx » secondary thrust coefficient parallel to primary nozzle axis 
cfy ■ secondary thrust coefficient normal to primary nozzle axis 

CMS B SECONDARY NOZZLE FLOY) COEFFICIENT 

CSTAkP"PR (MARY FEOTi CHARACTERISTIC VELOCITY (FEET PER SECOND) 

CSTAh ^.SECONDARY FLO* CHARACTERISTIC VELOCITY (FEET PER SECOND) 

DP ■ DIAMETER OF THE PRIMARY NOZZLE THROAT IFEET) 

EPS . SECONDARY INJECTION ANCLE (DECREES) 

CAMP" PRIMARY RATIO OF SPECIFIC HEATS 
CAMS" SECONDARY RATIO OF SPECIFIC HEATS 

CIMANC • G I MOAL ANCLE (RADIANS) 

PED • PRIMARY EXIT OIAMETER (FEET) 

PM a LOCAL PRIMARY FLOW MACH NUMBER 

PSL ■ PRIMARY NOZZLE LOCAL WALL STATIC PRESSURE (PSFA) 

PTP ■ PRIMARY TOTAL PRESSURE (PSFA) 

PTS « SECONDARY TOTAL PRESSURE (PSFA) 

S ■ INTERACTION DISTANCE IFEET) _ 

sh • secondary injection macm number 

SMI a INITIAL SECONDARY M A SSF L 0* R A T l 0 
sm i nc « secondary massflow increment 

SMMAX B MAXIMUM SECONDARY mASSFLOW 

SMR a SECONDARY TO PRIMARY MaSSFLO* RATE RATIO 

SPR B SECONDARY To PRIMARY INJECTION TOTAL PRESSURE RATIO 

T a SMK" AF 

THETaC b CIMBAL ANCLE (DEGREES) 

TIP a PRIMARY TOTAL TEmPERATUkE (DEGREES RaNK(NE) 

tts a seconoary total temperature (Degrees rankinej 

TU B bLAST WAVE ACCUMULATION OF TERMS 

VNP a VACUHN MOMENTUM SPECIFIC IMPULSE EFFICENCY (PRIMARY) 

V NS B VACUMN MOMENTUM SPECIFIC IMPULSE EFFJCENCY (SECONDARY) 

W B umeca 

*1 B BLAST WAVE OMEGA ONE DIMENSION LESS ENERGY 
W 2 B BLAST WAVE OMEGA TWO DIMENSION LESS ENERGY 
W3 b BLAST WAVE OMEGA THREE DIMENSIONLESS ENERGY 
WMP a PRIMARY GAS MOLECULAR WEIGHT 

wms a secondary gas molecular weight 

WPR 1 a PRIMARY MASSFLQW rate (POUNDS PER SECOND) 
w SEC ■ secondary MASSFEUW RATE (POUNDS PER SECOND) 

XL ■ RATIO OF DISTANCE FROM PRIMARY NOZZLE THROAT TO INJECTION POINT 
DIVIDED BY PRIMARY NOZZLE LENGTH 


REAL L 

READ(S,IOI SPH,XL,EPS,SM,CMS,VNS,PTP 
READ l *> , l 2 I ARATIO, GAMS, VNP, COP. SMI, PSL 
READ (S , I 9 ) WMP,*MS, T TP , WPR I 
RE AD ( b « 16) TTS, GAMP, PM, S. AS 


C-2 
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UO 1 HI 

6 M , 

Uu 1 *1 1 

6S* 

Ju 1 M 1 

_ 66* 

UOI M 1 

67* 

UOIMI 

68* 

uU 1 H 7 

69* 

U0I50 

7u* 

U0IS1 

7 1. 

UO l 52 

72. 

UU153 

73. 

ODIUM 

7 H * 

UOI 56 

75. 

UU 1 5 7 

76 . 

00IS7 

77. 

UO 1 6u 

7 8. 

UU 1 6 1 

7 V • 

UU 1 6/ 

80* 

UU 1 6/ 

81* 

UU 1 63 

82* 

UO 1 6 H 

»0 1 AG 

UU 1 6S 

83. 

0 0 1 6 6 

8 M • 

00167 

8 5* 

00 170 

86. 

0017 1 

87* 

uO l / I 

88* 

0017/ 

a?« 

UOI 73 

90* 

UOI 7M 

?1« 

UU 17 5 

92* 

00176 

93. 

UOI 77 

9 M . 

UO 2 0 0 | 

95. 

U020u 

96. 

00 2 0 1 

97. 

UU2U2 

98. 

00203 

99* 

002 0 H 

100* 

00205 

101 * 

00205 

102* 

UU 2 U 6 

103* 

00 2 0 7 

1 0M« 

00 2 1 0 

105. 

uO 2 1 1 

106* 

002 1 2 

107* 

002 1 2 

108. 

U02 1 2 

109. 

002 1 2 

110* 

002 1 3 

Ill* 

002 1 H 

112* 

Up 2 1 5 

113. 

UU 2 1 6 

1 IM» 

002 16 

1 IS* 

•302 1 6 

116. 

JO 2 l 6 

117* 

J021 7 

lie. 


READ ( 5 , I8ICFP, ALPHA #SM INC, SHm AX 

7 CFMP CALCULATION 

SMR ■ SMI - SM1NC 
[COUNT • 0 
IOU CONTINUE 

ALPHA - ALPHA/S7.29S82 
SHK « SMR « SMINC 
1FISMK ,GT, SMMAX) GO TO IQt 

CGAMS«(6AMS»»«S)»( (2.0/IGAMS.l . ) > • • ( ( G A M S * 1 , )/(2.«(GAMS-l. I I I ) 
CGAMP « ( (GAMP ) »..S).( (2.0/IG AMP ♦! ,0)1 •• I I GAMP* 1 .0)/ I 2 • 0* I GAMP- I • 
D) ) ) 

C S T ARP • SORT! ( 1 S M H • 0 / WMP ) • T T P I / C G AMP 
CSTAKO-SQRT ( ( I 5 M M . 0 / A MS > • T T S >/ C G A MS 

CFVO«(CGAMS)*( I .C.GAMS»SM«.2»0)/tSQRT(6AMS) )*(SM) .([..((GAMS-!.) 
l/2.0.Sll».2.0))...5 
ASECaWPR l .SMR 

t c • the test for eoualiTt between non-integers may not be meaningful. 

|F( <AS.NE.O»0> .and*. ( KOUNT. EO. OIL <*0 TO III 

PTS > SPR.pTP 

TERM 1 .SORT ( GAMS*«MS/ ISMM. 1 
TtRM2»PTS/SQRT|TTS) 

TERMJ-SM/I I .♦( (GAMS-1 . 1/2, >»SM»S«)»»I I G AMS* I • ) ✓ I 2 t • I G AMS- I . 
I))) 

AS ■ A5EC/ < TERM I »TERm2.TERm3.32.2> 

aseff»as«cms 

ICOUNf • 1C0UNT ♦ 1 

GO TO I 12 

111 ASEFF « A S • C MS 

TERM 1 "SORT (GAMS* AMS/ I SHI. I 

TERM3*SM/( I .♦( (GAMS-I . I /2.J_.SM*SM » I I G AMS* 1 . I / I 2 , • I G AMS- 1 . 

Ill) 

PTS ■ (MSEC* SORT ITTS) (/(ASEFF*TERM1»TErM3*32#2) 

112 CONTINUE 

TERM? - SURT ( GAMP*aMP/ )SMM, > 

TERMS - PTP/SQHT ( TTP I 

TERM* ■ 1,0/1 l .♦ < I GAMP- I . ) /2, ) • 1 , 1 •• I I GAMP*1 , 1/ l 2,* ( GAMP-I . 

tin ' 

AP ■ APR|/(TERMM»TERM5.TERM6*32,2) 

ASTARP ■ AP • COP 
DP ■ SORT < H. • A5TARP/3, 1N15?) 

0* ( PSL/PTP>*( 1 ./SMR ) *(CSTARP/CSTARsl*t ( ASEFF)/(CDP«ASTARP)) 
CFMP-tVNSI*(CFVO»-0 


CFMN CALCULATION 

A.(3.|HIS9)»(CGAMSI»((.S)-I EPS/ 180,0) I 
B-(CMS)*S0RT (GAMS)*SM 

C-SQRTt 1 .O+tGAMS-1 <01* ( SM* • 2 . 0 ) / ( 2 , 0 ) >_ 
CFMN* ( A/ ( B*C J ) - 0 / T A N ( £ P 5 / 5 7 » 2 9 5 8 2 ) " 

theory i induced forces (high smr> 

TU»( l . 17S) • ( GAMP »-l ,05 


C-3 


010 
01 1 


006 
012 
01? 
020 
02 M 

025 

026 
030 
0?0 

060 


065 

076 

077 

078 
07? 
080 
070 


08 1 
082 
083 
08? 

085 

086 
087 

1 20 
I 10 

I 30 

1 MO 
ISO 
160 
170 

180 

190 
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0021 7 

1 19* 


U02 1 / 

120* 


00 2 1 7 

121. 

C 

UO 2 2 u 

122. 


UO 2 2 1 

123. 


00222 

129. 


UU2 2 J 

125. 


00223 

126. 

C 

0022 J 

127. 

c 

00223 

128. 

c 

0022 9 

12V. 


0022b 

1 30* 


0 0 2 2 6 

I3|. 


00 2 2 6 

132. 


002 2 6 

133. 


00 2 2 6 

139. 

c 

00226 

13b. 

c 

0U226 

1 3o* 

c 

00227 

137. 


0022 7 

138. 

c 

0022/ 

139. 

c 

0022 / 

190. 

c 

00 2 3 o 

19 1* 


0 0 2 3 0 

192* 

c 

00230 

193. 

c 

00 2 3 u 

199. 

c 

0023 1 

1 9b. 


0023 1 

196. 

c 

0023 1 

19 7. 

c 

0023 1 

198* 

c 

00232 

199. 


00232 

1 bU. 


00232 , 

151. 


00232 

152. 

c 

00232 

153. 

c 

OU232 

159. 

c 

00233 

155. 


00239 

156* 


U023b 

15 7* 


00236 

158. 


0023 7 

159. 


0 0 2 9 0 

1 60. 


0 0 2 9 0 

161* 

c 

OU29J 

162* 

c 

00290 

163* 

c 

0029 1 

169. 


002 9 2 

1 6b. 


00 2 9.-. 

1 66 » 

c 

0 0 2 9 2 

1 67. 

c 

002 9 2 

168* 

c 

002 9 • 

169. 


JO 2 3 9 

170* 


u 0 2 9 5 

17 1* 


00 2 9 6 

1 72* 


0029/ 

173* 


J U 2 5 l- 

179. 



OMEGA CALCULATIONS, OMEGA ONE FIRST 

E» I I . O* I G AMP-T.0) * < I PM*»2. 0 1 /2.0IT ~ 

F«(CSTARP/CSTARS> 

G«( (S.0*GAMP>-3.0>/( ( 8,0»GAMp»-e.0) , 

» I • 1 ( P M • • l .25).(F».,75| I / C E • * 6 1 

OMEGA TWO CALCULATIONS ; _ ; 

H«( GAMP I « (GAMP-t ,U> • (PM»»2.U> 

P« ( GAMS/ ( GAMS-1 .01 I 

«2« * 5* ( I •U/H|*(P*E)/C6AMP»PM«*2«0I*I ( C G AMS/C G AMP I . * 2 . 0 ) • I I *0/F)*» 
I2.0«|Sm/Pm)*S 0KT(E|/C •IS«HT(GAMS/GAhP) ).(CG aMS/CGAMP).( I ,0/F )• 
2C0SIEPS/57. 29582) 

OnLuA CALCULATION 


W*W I • ( W 2 » » , 7 5 ) 

CF1NI CALCULATION 

~ . . ^ ~ 

C F INI" ( TUI • « Wl • I SQRTI s/ I SORT ( COP ) • ( OP) )| I • ( I .O/SMR) • ( F.. .25 ) 

THEOKV 2 INDUCED FOHCES (LOW SMR) 

« - S«RT|6aMPI*CGAMP*(GAMP-I,0»*(Pm*0.0) 

OMEGA THNEt CALCULATION 

W3-W2 ♦( 1 .0/ ( (GAMP- I .0) •PM«*2.0> I • ( WMP/WMSl • ( I .-( I GAMP- l .»/ I 6 AMS- l 
1,M*GA MS/GaMP) ♦ (GAMP - l.UI/GAMP ♦ ( SM/PM ) • ( SORT I E ) /C I .SORT I 
2GAMS/GAMP)»(CGAMS/CGAMPI*(|.0/FI»C0S|EPS/&7.29S82) 

CFIN2 CALCULATION 

CF1N2 •t«»F»V,3l/M (Pm»*2.0I-I.0»»«E> >••«& 

PED • S6IHT(ARAT10«DP»»2.0I 

L ■ 1.0 - «6S«(S/PE0)»*1 .25 

CFIN1».7*L*CFIN1 

Cf I N2 ■ l t 0* L *C F I N 2 

CFIN -aM1N( (CF1N1 ,CF I N 2 » 

AMPLIFICATION FACT0« C ALCUL AT I ON I I NTER ACT I ON CONTRIBUTION FIRSTI 

' CFIX ■ CFI(.»SIN|ALPHaJ - 

CF I Y » CF IN»C0S(ALPHA) 

AMPLIFICATION factor (MOMENTUM CONTRIBUTION) 

H « alpha ♦ EPS/S7. 29582 | 

C F M X ■ -CFMP*COS I R I ♦ CFMN.SlNIR) “ ~~ 

C F M Y ■ CFmP.SIN(R) ♦ CFMN.COStRI 

CF X ■ CFMX ♦ CFIX 

CFY • CFMY ♦ CF I Y 

AF • (CFY«wSEC*CSTaRS)«(I.O/SMR)/ICFP**PRI*CSTARP ♦ CFX*WSEC» 


C-4 


200 

‘ 210 
220 
23U 
290 

2S0 

260 

270 

280 

290 

300 

310 

320 

330 

390 

350 

360 

370 

380 

390 

900 

910 

920 
922 
990 
~ 950 
960 
930 

970 

980 

990 

500 

520 
530 
590 
S 50 
560 
570 
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0 U 2 5 0 

1 7S* 

002b 1 

176* 

0(1252 

1 7 7* 

00253 

1 78* 

002 5 9 

1 79* 

002 b S 

1 80* 

00259 

181* 

0025b 

. 162* 

0l>25 7 

183* 

00260 

las* 

0026 7 

185* 

00270 

1 86* 

00272 

187* 

00273 

188* 

00302 

189* 

00301 

1 9U* 

O0305 

191* 

003 1 S 

192* 

00316 

193* 

0031 7 

19S* 

003 1 7 

195* 

003 1 7 

196* 

00322 

197* 

00323 

198* 

00 3 2 9 

199* 

00325 

2 06* 

00 3 2 6 

201* 

003 2/ 

2 U 2 • 

00 3 3 0 

203* 

0003 1 

20S • 

00332 

205* 


END OF 


c 

c 


t 

c 


ICSTARSl 

T«SMR*AF 

GIMAnG ■ ATAN(T ) 

the tag • 6i hang • s7 , 29582 

ALPHA • ALPHA*S7. 29582 


580 

581 

582 

583 
58S 


WRITE (6. 200) - 

2 «Q FORMAT ( (HI . , iy*,2HAS. 12Ki.lHCFHP.II>. HHCFMN t Hx t 6HTHETA6.yX.3HSMR) , 

WRITE (6, 201 lAS.CFMP.CFMN. THETAG iSMR 

201 FORMAT ( 1H0.9X , IPSE 15.5) 

WRITEU.202) 

202 FORMAT! 1 HO , 1 8 X , 5HC F I N 1 . 1 OX , 5 HCF I N2 , 1 OX , S HCF I N , 1 2X , 2H AF , I OX , 3HPTS ) 

WR 1 TE ( 6 » 20 1 ) CFIN1 , CF I N2 i CF 1 N , AF , P T S 

WR 1 TE ( 6 t 203 ) 

203 FORMA T ( I HO . 1 8X . 3HCF* , I3X ,3hCF T , I I X , SHCFMXTifx .SHCFMY ,T6* .5HALPHAT “ 
WR I TE ( 6 • 20 1 ) CFX.CFY.CFMX.CFMV, ALPHA 

WRITEU.20SI 

20S FORMAT! 1H0, 1«X, SHFS/FA) 

WR 1 TE < 6 • 20 1 IT 


10 FORMAT ( 7E1U«5) 
12 FORMAT ( 7 1 1 0 • 5 1 
IS FORMa T ( 6E 1 u • 5 ) 
U FORMA T ( 7E 10*5 ) 
18 FORMAT (5E1Q. 5) 
GO TO 100 
101 CONTINUE 
STOP 
END 


590 

800 


COMPILATIONS 


I diagnostics. 
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Appendix D 

NON-FLOW THROUGH PEBBLE 
BED HEATER PERFORMANCE 


P ' 1 
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Appendix D 


The performance of the non-flow through pebble bed heater described 
in Section 4.3 of this report is calculated by applying the first law of thermo- 
dynamics and the perfect gas equation of state. Applying the first law, one 
obtains: 


The decrease of stored 


The stored energy of 


The flow work of 

energy of CF^ in the tank 


the gas leaving the 
tank 

T 

the gas leaving 
the tank 





- 


The above equation written in differential form is 


mdu = p Pdm 


The specific internal energy of the CF^ may be expressed as 


u = .9313 pp 


The energy balance becomes 


= 0.9313 -^EXl 
m pv 


or 


£n 


m. 


m 


u 


0.9313 in 


Pf v f 

Pi Vi 
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where the subscripts f and i refer to the final and initial conditions in the 
heater, respectively. Further manipulation yields 



This equation is somewhat conservative as the amount of heat energy added 
to the CF^ through heat transfer from the pebbles was neglected. 


A sample calculation using the above equation and the equation of state 
follows. 


C ons ide r: 


and 


p^ = 500 psia; p^ = 300 psia; 


= 20 pounds of CF^ 



15.64 pounds 


Now for a set flow rate of say 2.5 pounds per second one would obtain a test 
"blow" of 2.5 seconds duration. Also since pv =RT and V = v/m. , one can 
solve for the required volume of the heater which for this case is 5.3 cubic 
feet. 

It is emphasized that the preceding analysis is restricted to 100% CF. 

o ^ 

heated to 600 F. Analogous results may be obtained for other gases and other 

test conditions. However, the appropriate thermodynamic constants applying 

to the particular gas or gas mixture being considered must be used. 
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